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&<^b¥aw&Qj  a         ABSTRACT 

This  thesis  presents  a  method  for  estimating  the  distribution  of 
stress  components  around  the  outer  circumference  of  a  cross  section  of 
a  pipe,  from  the  readings  of  strain  gage  elements  arbitrarily  positioned 
and  oriented  around  this  circumference.   A  least-squares  procedure  is 
used  to  obtain  best  estimates  of  the  coefficients  of  Fourier  expansions 
describing  such  distributions.  A  digital  computer  program  was  developed 
for  applying  the  method.   Data  for  testing  the  method  and  program  were 
generated  by  a  computer  program  using  the  best  available  theory  of  stress 
analysis  in  pipes.  Methods  for  adding  random  errors  to  the  data  were 
adapted  and  used  for  closer  simulation  of  actual  situations. 

A  second  problem  treated  in  this  thesis  is  that  of  inferring  the 
loading  acting  through  the  cross  section  of  a  straight  pipe  of  concen- 
tric bore  from  known  stresses  at  points  on  the  external  surface  of  the 
cross  section  which  is  presumed  to  be  distant  from  stress  concentra- 
tions.  It  is  shown  that  this  inference  can  be  made  from  stresses  at 
only  three  points  of  the  cross  section.  A  digital  computer  program  was 
developed  to  do  this. 
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1.    Introduction. 

In  a  thesis  presented  to  the  Naval  Postgraduate  School  in  May,  1966, 
LT  Joseph  W.  Koch,  Jr.,  USN,  did  pioneer  work  on  the  problem  of  obtain- 
ing optimal  inferences  of  stresses  in  submarine  main  sea  water  piping, 
from  strain  gage  readings  taken  during  dockside  hydrostatic  tests  and 
during  submerged  operations.   This  work  provided  a  sequential  procedure 
for  determining  such  stresses,  but  it  also  indicated  the  need  for  many 
related  investigations. 

Two  of  these  are  treated  in  the  present  thesis.   The  first  is  that 
of  inferring,  in  an  optimal  way,  the  distributions  of  stresses  around 
the  outer  circumference  of  a  given  cross  section  of  pipe  from  observa- 
tions of  strain  obtained  from  an  arbitrary  distribution  of  strain  gage 
elements  around  that  circumference.   There  are  several  reasons  for  being 
interested  in  this  problem.   One  reason  is  that  even  though  the  gage 
elements  may  be  arranged  in  special  patterns  intended  to  simplify  deter- 
mining stresses  (for  example  three-gage  rosette  arrangements),  failure 
of  any  of  the  elements  might  make  it  difficult  to  make  effective  use  of 
valid  information  obtained  from  closely  related  gages.   Another  reason 
is  the  ability  to  infer  stresses  from  arbitrarily  arranged  gages,  which 
is  a  logical  step  in  investigating  the  most  effective  and  most  economi- 
cal arrangements  of  gage  elements. 

Briefly,  this  problem  is  solved  herein  by  constructing  Fourier 
series  representing  the  desired  stress  components.   Theoretical  consider- 
ations are  used  to  evaluate  the  coefficients  of  the  series  in  an  optimal 
way,  minimizing  the  sum  of  the  squares  of  the  residuals.   A  digital  com- 
puter program  was  developed  to  support  the  theory. 


The  second  problem  treated  here  is  that  of  inferring  the  force  and 
moment  components  applied  through  a  cross  section  of  pipe  and  the  inter- 
nal pressure  acting  inside  the  pipe.   It  is  shown  that  these  values  may 
be  inferred  from  a  knowledge  of  the  state  of  stress  at  only  three  abritrar- 
ily  chosen  points  at  the  external  circumference  of  the  cross  section,  and 
a  digital  computer  program  for  making  the  corresponding  evaluations  was 
developed. 

The  solution  of  this  problem  assumes  that  the  stress  distribution  is 
that  given  by  what  may  be  regarded  as  the  best  available  current  theory 
for  the  state  of  stress  in  a  section  of  uniform  pipe  having  concentric 
bore  and  remote  from  stress  concentrations  and  other  disturbances  in 
the  stress  pattern. 

A  digital  computer  program  to  generate  numerical  data  from  the  best 
available  theory  was  developed.  This  data  was  used  to  test  and  verify 
the  theory  and  the  computer  programs.   Suggestions  of  how  to  extend  the 
present  work  are  given  in  Section  4. 

The  author  would  like  to  express  his  sincere  appreciation  to  Dr. 
John  E.  Brock,  Professor  of  Mechanical  Engineering,  Naval  Postgraduate 
School,  for  his  continued  patience,  efforts  and  most  capable  guidance 
while  acting  as  faculty  advisor  in  this  work. 


2.    Least-Squares  Estimation  of  Stress  Distributions  at  the  external 
Circumference  of  a  Pipe  by  Means  of  Fourier  Expansions,  Using 
Strain  Gage  Data. 

The  usual  method  of  estimating  the  stress  distribution  at 
the  surface  of  a  body  is  by  taking  readings  of  strain  gages  grouped  in 
rosettes,  to  determine  local  principal  stresses.   Unless  the  strain 
rosettes  have  a  redundance  of  gage  elements,  when  one  or  more  of  them 
fail,  it  is  difficult  to  make  effective  use  of  the  information  provided 
by  the  elements  which  do  not  fail. 

The  work  that  follows  presents  a  method  to  estimate  the  stress 
distribution  around  an  external  circumference  of  a  cross  section  of  pipe, 
along  a  sequence  of  strain  gage  elements,  which  may  or  may  not  be  group- 
ed in  rosettes.   If  one  or  more  gage  elements  fail,  the  information 
from  the  others  is  still  useful. 

It  is  assumed  that  strain  gage  elements  are  arbitrarily  placed 
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Fig.  2.1. 
and  arbitrarly  oriented  around  the  external  circumference  of  a  cross 
section  of  pipe  so  as  to  measure  the  strain  produced  by  application  of 


loading.   Assume  that  there  are  N  such  elements  as  indicated  in  Fig. 
2.1. 

Define  two  orthogonal  directions  at  the  cross  section,  one  normal 
to  it,  called  longitudinal  and  represented  by  the  letter  L,  and  the 
other  tangential  to  the  external  circumference,  called  tangential  and 
represented  by  T;  Fig.  2.1.b.   Also  define  the  angle  9  as  the  angle  be- 
tween the  vertical  line  that  passes  downward  through  the  center  of  the 
cross  section,  and  the  radial  line  that  passes  through  any  specified 
point  of  the  external  circumference,  measured  as  indicated  in  Fig.  2.1.b. 

Measure  the  angular  orientation  of  each  gage  element  in  the  counter 
clockwise  direction  from  the  tangential  axis  T,  and  call  this  angle  0. 

A  set  of  parameter-pairs  can  then  be  formed 


(9.   ,  0.) 


i  =  1,  2,  ...  ,  N 


which  describe  the  position  and  orientation  of  each  element,  Fig.  2.2, 


C\       TeT  e>~v  o  7T        o\        360° 

1  »    •    2  l  x   N-l  i        N\ 


Fig.  2.2. 
At  each  one  of  the  gage  locations,  there  are  tangential,  longi 
tudinal,  and  shear  strains,  (€  )  ,  (C  )., 

XX         u  1 

and  (  €.   _,.).  respectively.   The  strains 

1  Li   1 

€   and  €   along  axis  y  and  z,  rotated 
y       z 


the  angle  0  from  the  tangential  direction 
Fig.  2.3.,  are  related  to  the  local 


Fig.  2.3. 
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strains  €  ,  €   and  €    by  the  strain  tensor  in  the  following  way 

1        Ij  1  Li 


y    yz 


yz    z 


cos  0    sin 


•sin  0    cos  0 


T    *  TL 


t  \ 


TL 


cos  0    -sin 


sin  0     cos  0 


The  corresponding  relations  for  the  strains  (  €  )   and  (  «  ) 
along  and  transverse  to  the  ith  strain  gage  axis,  are  then 

(%}i   "   (*T)iSinX  +  <*L>ic08\  +  <€TL)i,in  0iC°8  0i    ^2.1.1) 

(«2)i   -   («T)1co»201  +  (6L)iSin20  -  («TL)18intf1cos  0t  (2.1.2) 

At  the  external  surface,  strains  are  related  to  the  corresponding 
stresses  as  follows 

*T  «  <sT  -  v  sL)  /E  (2.2.1) 

«  ,   -  (sT  -  v  sj  /E  (2.2.2) 


TL  -  STL  '  G 


(2.2.3) 


where   s  is  stress 

v  is  Poisson's  ratio 

E  is  Young's  modulus  of  elasticity 

G  is  the  shear  modulus  of  elasticity 
Substituting  these  expressions  for  the  strains  in  the  relations 
(2.1)  for  the  strains  along  and  transverse  to  the  ith  gage  element  and 
simplifying 

(€    ).    -   <*T>i11   "    (v  +  Dsin^J/E  +  (sL)i[l   -    (v  +  Dcos^J/E 


+(sTL).sin  2  0,7(20) 


(2.3.1) 


(€z)i  =   (sT)[l   -    (v  +  Dcos^J/E  +  (bl)±[1   -    (v  +  l)Bin\]/E 

-    (sTL).sin2  0./(2G)  (2.3.2) 

A  strain  gage  element  is  sensitive  not  only  to  the  strain  along 
its  own  axis  but  also,  to  a  small  extent,  to  the  strain  at  right  angles 
to  this  axis.   The  transverse  sensitivity  is  usually  hard  to  obtain.   It 
is  rarely  available  from  the  manufacturer  and  may  have  to  be  inferred  from 
an  experiment.   Nevertheless,  it  may  be  of  importance  and  it  is  included 
in  the  theoretical  development  that  follows. 

For  a  particular  gage  element  properly  mounted  on  a  given  definite 
material,  the  relative  change  in  electrical  resistance,  which  is  sensed 
by  the  strain  indicator  (or  recorder) ,  is  given  by  a  formula  of  the  sort 

Ar  -  F   €     +   F  €  (2.4) 

—    y  y     z     z 

where  R  is  the  electrical  resistance  of  the  element,  6   and  £  are 

y       z 

the  strains  along  and  perpendicular  to  its  axis,  and  F  and  F  are  co- 
efficients which  depend  on  the  characteristics  of  the  gage  and  the  pro- 
perties of  the  material  upon  which  it  is  mounted. 
The  ratio 

k  =  F  /  F  (2.5) 

z    y 

is  called  the  transverse  sensitivity  (factor)  of  the  gage,  and  the 

quantity 

F1  =  (1  -  k2  )  F 

y 

is  called  the  gage  factor  (which  should  be  set  on  the  sensor  instrument). 
F'  is  the  ratio  of  unit  change  in  electric  resistance  to  the  quantity  e, 
(indicated  strain)  sensed  by  the  instrumentation. 
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Thus  e   =      (     ^     )    /  F1 

K 

-   (  F      6       +  F      €     )/[(l   -  k2)F   ] 
yy         z       z'  y 

=   (   €y  +     k€z   )/(l   -  k2)  (2.7) 

This   is  what  the  instrumentation  should   indicate   if   there  were  no 
errors.      Clearly  for  a  gage  with  no  transverse  sensitivity,   k  ■     0,   and 
e  *    € 

y 

After  substituting  the  expressions  (2.3)  for  (€  ).  and  (  €•  ) .  in 
equation  (2.7)  and  simplifying,  the  theoretical  indicated  strain  e 
sensed  along  the  ith  strain  gage  element  is  given  by 

e.   =  (sT)1[(l  -  kv)  -  (1  -  k)(l  +  v)sin20i]/[E(l  -  k2)]  + 

+  (sL)i[(l  -  kv)  -  (1  -  k)(l  +  v)cos20i]/[E(l  -  k2)2  + 

+  <«TL>1<1   "  k)sin2(Z>i/[2G(l   -  k2)]  (2.8) 

Let  the  coefficients  of  the  stresses  in  the  above  expression  be 
represented  by 

AC^)  =  A±   =  [(1  -  kv)  -  (1  -  k)(l  4v)sin20.]/[E(l  -  k2)]   (2.9.1) 

B(0t)   =  B.    -   [(1   -  kv)    -    (1   -  k)(l  +  v)cos20i]/[E(l   -  k2)] 

(2.9.2) 

C(0±)   -  C±  =  (1   -  k)sin20i/[2G(l  -  k2)]  (2.9.3) 

Then  equation   (2.8)   can  be  put  in  the   form 

Any  one  of  the  stresses  (s_,) .  ,  (s_)  and  (sa<T)  has  continuous  and 
finite  values  around  the  external  circumference  of  the  pipe  cross 
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stress 


Fig.  2.4. 
section,  and  can  be  represented  as  in  Fig.  2.4. 

The  functional  relation  between  stress  and  the  angular  position 
0  can  be  expanded  in  an  infinite  Fourier  series 

(s__)  sin  n9. 


<Vi  - 


<8L>i  " 


(,TL>1  * 


where 


^~*   (smA)  cos  n6J  + 
•*-o   v  TA'n      i 


H;  o 


TB  n 


(s_.)  cos  n0J  + 
LA  n      i 


S  ^TlA"8  nei  + 


^-  (s_„)  sin  n9. 
"=«  v  LB'n       i 


f-]   (Wnaln  n9i 


(2.11.1) 
(2.11.2) 
(2.11.3) 


(Vn  '  ^n  *   <sIA>n   »   (8LB>n   •   (sSHA>n   '   and  <8SHB)n  are 
as  yet  unknown  numerical  coefficients. 

Sufficiently  accurate  values  can  be  obtained  if  the  expansion  is 
carried  to  a  finite  number  of  terms,  say  P  terms,  in  which  case  the 
index  n  would  reach  a  maximum  value  of  P/2.   Call  this  value  M.  Then 

(s,D)_sin  n0,      (2.12.1) 


<Vt   » 


n-o 


(sm.)  cos  n9.  + 
TA'n      i 


Is  a 


TB'n 


The  conventional  expansion  of  a  Fourier  series  is  a  constant  term 
followed  by  a  summation  from  1  to  co  .   The  same  is  obtained  by  eliminat- 
ing the  constant  term  and  making  summations  from  0  to  ©o  .   This  latter 
form  will  be  used  in  this  work  to  facilitate  mathematical  manipulations 
and  later  programming  for  a  digital  computer. 
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M  H 

(st>-  CZ>  *£      (sia)  cos  n0„  +    2*   (s__)  sin  0.       (2.12.2) 

L  1  ~~  mbo    LA  n      l      Hsc    LB  n     i  ' 

/*  M 

<8TL>i  «     SI  (sTLA)nCos  n0.  +    C   («TLB)nsin  n^     (2.12.3) 

Note  that  this  expansion  contains  6(M  +  1)  Fourier  coefficients. 
If  these  expressions  for  the  stresses  are  substituted  into  equation 
(2.10)  for  the  strain  e  ,  theoretically  sensed  by  the  ith  strain  gage 


element, 


M  M 


e,  C2        A.  [    -Ei    (s  )  cos  n0  +  S  ...']  + 
l  i    >,-0  TA  n      l    *,g0 

+  B.[   2ES    ...   +       JC        ...    ]   +        (2.13) 

+  C[  aS  ...+   2X   (s   )  sin  n0  ] 
i   x^c         *-0    TLB  n      i 

Call  e*  the  strain  actually  measured  by  the  ith  strain  gage  ele- 
ment, and  E.  the  difference  between  that  measured  strain  and  e.,  the 
l  l 

strain  that  theoretically  it  should  have  measured,  so  that  this  error 

E.  is  expressed  as 
i 

E.  =  e.  -  e*  (2.14) 

ill 

Let  R  be  the  sum  of  the  squares  of  the  errors,  that  is 


R  = 

Z2 

(E,)2 

A/ 

*B* 

\   ^/ 

- 

£. 

A.CI 

ci: 

e*)   2 

<V2 

-     2 

SC       e.e*  + 
*.*#         i  i 

,?> 

(2.15) 

To  minimize  R,  partial  derivatives  respect  to  each  of  the  6(M  +  1) 
Fourier  coefficients  are  taken  and  set  equal  to  zero: 

~%  R    ■  0 
where  Dk  =  (sTA)k  ,  (s^  ,  (s^  ,  (8LB)k  ,  (s^  ,  (s^  in  turn 
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for  k  -  0,  1,  2,  . ..,  M 

Taking  derivatives  of  R  from  expression  (2.15) 


4  R        '0   =  2     5£ 

e. 

i 

^ 

"b 

e. 

i 

Dk 

-   2 

N 

->Dk 

Simplifying  this 

A                "*    e. 
/■I       *           ^    Dk 

V 

el 

*Ek 

*  ^tA 


^6i 


^  (STlA 

^6i 


(2.16) 


(2.17) 


Taking  derivatives  of  e  from  expression  (2.13) 

^  e 

i -  Aicos  k0i   for  k  -  0,  1,  2,  ...  ,  M        (2.18.1) 


^~yk  ,   A1.ink91  (2.18.2) 

~»6i 
^Ti-?k   -  B^oskO.  (2.18.3) 

2)  ei 


1  (*l£>k       =  Bisin  k9i  (2.18.4) 

*  ei 

=  C^cos  ^  (2.18.5) 


-*(sTLB)k  -  C.sin  k0i  (2.18.6) 

Note  that  there  are  6(M  +  1)  derivatives  of  this  form.   Substitut- 
ing in  equation  (2.17)  the  expression  for  e.  from  equation  (2.13) 

friAi  L  ?E  (Wnsln n9i +  •••J+  ••• +  ci£  ••• 

•  •  •  +     f-     <  W„sin  n9iJr  JdT     "    S   «f   -^ 

*SO  -JJ     "     k  4,'-,  w       k 

(2.19) 
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which   is   the   same  as 


c  (i.)  r£  a. erne.  i&  7  +  zjcsla  [sa^ -»*.?&  i 


for  k  -  0,  1,  2,  ...  ,  M 

When  the  index  n  in  the  Fourier  expansions  is  equal  to  zero,  the 

terms  containing  sin(n0.)  are  zero.  This  makes  the  coefficients  (s__)-., 

1  Id  (J 

(s..,)-.  and  (s--^)-  also  zero,  reducing  the  Fourier  coefficients  to  (6M  + 
LB  u        TLB  0  ^ 

o  e . 

3)  and  the  same  number  of  derivatives  of  the  form   ^   .  There  will 

k 
be  (6M  +  3)  equations  in  *each  of  which  appear  all  the  Fourier  coeffici- 
ents. This  makes  a  set  of  (6M  +  3)  equations  with  the  same  number  of 
unknowns,  (the  Fourier  coefficients),  which  can  be  solved. 

Once  the  Fourier  coefficients  are  solved  for,  they  can  be  substitut- 
ed in  expressions  (2.7)  for  the  stresses.   The  stress  distributions  s_, 
8T  and  s   are  defined  then  as  functions  of  the  angle  9  over  the  com- 

L  Li-i 

plete  external  circumference  of  the  cross  section. 

A  digital  computer  program  was  developed  to  apply  the  theory  of  the 
preceding  text.   It  was  tested  using  a  mathematical  model  of  a  canti- 
levered  pipe.   The  program  is  called  TERESITA  and  is  presented  in  Ap- 
pendix B  of  this  thesis. 

It  was  assumed  that  strain  gages  having  a  zero  transverse  sensiti- 
vity factor  were  placed  at  the  external  circumference  of  cross  section 
No.  1,  see  Fig.  A.l.  in  Appendix  A.  Four  arrangements  of  gage  elements 
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and  several  numbers  of  Fourier  coefficients  were  used. 

The  first  arrangement  consisted  of  twelve  strain  gage  rosettes 
equally  spaced  30° from  each  other  around  the  cross  section.   Each  rosette 
has  three  gage  elements  oriented  at  0?  45°and  90° from  the  positive  tangen- 
tial axis.   Runs  were  made  with  11,  9,  7,  5  and  3  Fourier  coefficients 
for  each  stress  component. 

The  second  arrangement  consisted  of  twelve  rosettes  spaced  as  above 
but  with  elements  oriented  at  0°,  60°,  and  120°.  Runs  were  made  for  the 
same  number  of  Fourier  coefficients. 

The  third  arrangement  consisted  of  eight  strain  rosettes  equally 
spaced  45°  from  each  other.   The  rosettes  had  four  gage  elements  orient- 
ed at  0°,  45°,  90°  and  135°  from  the  positive  tangential  axis.   Runs  were 
made  for  the  same  number  of  Fourier  coefficients. 

As  input  data  for  testing  the  program,  stress  distributions  pro- 
duced by  the  loading  were  evaluated  at  the  external  circumference  of 
cross  section  No.  1.   From  the  stresses,  the  theoretical  strain  that  the 
gage  element  would  sense,  was  computed.   The  stress  distributions  re- 
covered by  the  program  agreed  in  all  cases  with  the  input  information. 
Results  using  eleven  Fourier  coefficients  (for  each  stress  component  s_, 
sT  and  s  )  are  almost  identical  with  those  using  only  three  coefficients. 

Additional  runs  were  made  using  the  second  arrangement  of  gages, 
but  eliminating  arbitrarily  the  readings  of  some  of  them.   This  is  what 
would  have  to  be  done  in  an  actual  case  if  some  gage  elements  were  known 
to  have  failed.   Results  are  satisfactorily  close  to  the  preceding  ones. 
Other  runs  were  made  using  a  strain  gage  transverse  sensitivity  factor 
of  1%.   The  program  recovered  likewise  the  input  information. 

The  fourth  gage  arrangement  had  one  gage  element  every  10°  around 
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the  circumference  with  different  orientation  for  each  gage.  Runs  were 
made  using  the  same  number  of  coefficients. 

All  the  runs  described  above  were  repeated  introducing  in  the  gage 
element  readings  small  random  errors  so  to  simulate  actual  situations. 
Calculated  loadings  differed  from  the  input  loadings  to  a  degree  which 
depended  in  each  case  on  the  amount  of  error  introduced. 

The  program  and  some  sample  numerical  results  of  the  runs  made  are 
presented  in  Appendix  B  of  this  thesis. 

A  more  detailed  discussion  of  the  effects  of  changing  the  number 
of  Fourier  coefficients  and  gage  element  arrangement  is  presented  in 
Appendix  E. 
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3.    Inference  of  the  Loading  Acting  at  a  Cross  Section  of  a  Pipe  of 
Concentric  Bore  from  Stresses  at  the  External  Circumference  of 
that  Cross  Section. 

If  stresses  are  known  at  a  sufficient  number  of  points  around 
the  external  circumference  of  the  cross  section  of  a  pipe  of  concentric 
bore  at  a  distance  from  stress  risers  and  concentrated  loads,  the  load- 
ing acting  at  the  cross  section  can  be  determined  from  theory. 

It  is  intended  to  analyse  the  effects  that  such  loads  produce 
at  a  general  point  P  under  the  conditions  specified  above.   From  this 
analysis,  general  conclusions  will  be  reached  which  lead  to  a  mathe- 
matical procedure  to  determine  the  loads. 

Consider  a  cross  section  of  pipe  as  specified  above,  Fig.  3.1. 
Choose  reference  axes  x,  y  and  z  with  the  center  point  0  as  origin.   Let 


l~qD^  sl 


I 


AA  view 


\    8TL 


(a)  (b) 

Fig.    3.1. 

i,  J  and  1c  be  unit  vectors  along  the  positive  x,  y  and  z  axis  so  to  form 

a  right  handed  triad.   Let  e"  be  a  unit  vector  along  the  outwards  direc- 

tion  of  the  radius  that  passes  through  the  point  P,  and  e   another  unit 

vector  so  that  e  ,  e   and  k  form  a  right  handed  triad.   The  angle  9  is 
R    T 
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measured  from  i  counterclockwise  to  e  .   Consider  the  differential  area 

R 

surrounding  the  point  P.  Assume  that  tensile  stresses  s  and  s   (along 
the  directions  of  e  and  k"  respectively)  and  shear  stress  s   ,  are  pro- 

■L  LLi 

duced  by  internal  pressure  p  and  by  a  general  force  F  and  a  general  mom- 
ent M  acting  on  the  cross  section  at  the  center  point  0. 

Let  the  components  of  F  and  M  in  the  e  ,  e  and  k  triad  be  (F  ,  F  , 

R    T  R    T 

F  )  and  (ML,  M_,  M  )  respectively,  and  in  the  i,  j,  k  triad,  (F  ,  F  ,  F  ) 
z       R   x   z  x   y   z 

and  (M  ,  M  ,  M  ). 
x   y   z 

From  strength  of  materials  considerations  reviewed  in  reference  5, 
the  stresses  produced  by  p,  F  and  M  at  the  external  point  P  are: 


sT=2pa     /(b     -a) 


,L  =   -^  b   /  I  +  Fz   / 


sTL  -  Mz  b   /    (2   I)   +  D  FT 


(3.1) 
(3.2) 
(3.3) 


where  a  *  internal  radius  of  the  cross  section 

b  ■  external  radius  of  the  cross  section 

A  =  cross  sectional  area 

I  =  moment  of  inertia  of  the  cross  sectional 
area  about  a  diameter 

D  =  [2  n  a2  +  (2n  -  2)b2]  /  [4I(  1  +  v) ] 

v  ■  Poisson's  ratio 

n  =  v  +  3/2 

Making  use  of  the  two  dimensional  coordinate  transformation  tensor 


cos  6     sin  0 


-sin  9 


cos  0 


M 


M 


M  can  be  expressed  as 
T 


M„,  -  -  M  sin  0  +  M  cos  9 
T      x         y 


(3.4) 
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Similarly, 


F_  =  -F  sin  9  +  F  cos  0  (3.5) 

T     x         y  v   ' 


The  expression  (3.2)  for  s  can  be  transformed  to 

Li 

sT  =  (1/A)F   +  (b/l)sin  0  M  -  (b/I)cos  9  M        (3.6) 
L        z  x  y 

Multiplying  through  by  (I/b) 

(l/b)s_   =  (l/(Ab))F   +  (sin  9)M  -  (cos  0)M       (3.7) 

Li  z  x  y 

Letting  (I/b)   =  g,  and  (l/(Ab))  =  h,  both  of  which  are  cross 
sectional  constants,  substituting  above  and  changing  sides 

h  F  +  sin  0  M   -  cos  0  M  -  g  s_  (3.8) 

z      x       y    L 

The  preceding  theory  is  applicable  to  a  point  on  the  exterior 
surface  of  a  piece  of  straight,  uniform,  concentric  pipe  at  a  suf- 
ficiently great  distance  from  its  ends  or  other  connections  which  might 
introduce  local  disturbances  into  the  stress  distribution.   Presuming 
that  these  conditions  are  satisfied,  we  next  turn  to  the  matter  of  try- 
ing to  infer  the  loading  conditions  at  the  cross  section,  i.e.,  F  ,  F  , 

x   y 

F  ,  M  ,  M  ,  M  and  p,  from  a  knowledge  of  the  stresses  at  several  points 
on  the  external  surface  of  the  cross  section. 

For  a  general  point  P,  8_  and  0  are  known  and  the  only  unknowns  in 
equation  (3.8)  are  F  ,  M  and  M  .   Taking  three  such  points  at  the  cross 
section,  the  three  unknowns  can  be  determined  proceeding  as  follows. 

Assume  that  the  angular  positions  for  the  points  are  0- ,  0.  and 
0-  and  the  axial  stresses  are  s. ,  s_  and  s_  respectively.  A  set  of 
equations  can  then  be  formed 

h  F  +s0.  M  -c9.  M  »  gs.  (3.9.1) 

z     1  x     1  y   °  1 

h  F     +  s0„  M     -  c0o  M     -  gs0  (3.9.2) 

z2x2y2 
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h  F  +sO,  M   -  c9.  M  -  gs. 
z     3  x     3  y   °  3 


(3.9.3) 


where  the  letters  s  and  c  with  no  subscripts  are  abreviations  for  sin 
and  cos  respectively.   These  abbreviations  will  be  used  throughout  this 
section. 

The  set  of  equations  (3.9)  is  solved  applying  Cramer's  rule 


Denominator  = 


s9, 


s9. 


s9. 


•c9. 


-c9. 


•c9. 


-h(  s(G1  -  92)  +  s(92  -  03)  +  s(03  -  0X))  (3.10) 


Numerator  of  F. 


gs. 


gs, 


gs. 


se. 


s9, 


-c9, 


-c9. 


s93    -c93 


=  -g(  Sls(92  -  93)  +  s2s(93  -  9X)  +  s3s(91  -  op) 

(3.11) 


Numerator  of  M 


h 

gSl 

-c01 

h 

gS2 

-c92 

h 

g83 

-c93 

-hg(  c91(s3  -  s2)  +  c92(Sl  -  s3)  +  c93(s2  -  s1)) 


(3.12) 


Numerator  of  M 


h 
h 

h 


s9. 


89, 


s9. 


gs. 


gs. 


gs. 


-hg  (s91(s3  -  s2)  +  s92(s1  -  s3)  +  s93(s2  -  s^) 


(3.13) 
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The  values  of  F  ,  M  and  M  are  found  by  dividing  the  respective 
numerator  by  the  denominator. 

The  remaining  components  of  the  vectors  F  and  M  appear  in  the 
expression  (3.3)  for  s__,  which  by  using  equation  (3.5)  can  be  trans- 
formed to 

s__   =   (b/2I)M   -  Dsin  9  F  +  Dcos  9  F  (3.14) 

TLz  x  y  N/ 

Multiplying  through  by  (1/D),  letting  (-1/D)  -  g*,  (-b/2ID)  =  h*, 

both  cross  sectional  constants,  substituting  in  (3.14)  and  changing 

sides 

h*  M   +  sin  9  F   -  cos  9  F   =  g*s„,T  (3.15) 

z  x  y    °     TL 

This  expression  is  of  the  same  form  as  (3.8)  which  led  to  the 

evaluation  of  F  ,  M  and  M  by  using  information  from  three  points.   In 

this  case  the  unknowns  are  M  ,  F   and  F  which  can  be  evaluated  by 

z       x      y 

using  the  same  three  points  and  applying  the  same  method.   In  this  form 
the  loads  F  and  M  are  completely  determined. 

Since  s   is  produced  by  the  internal  pressure  p  alone,  it  can  be 
solved  for  from  equation  (3.1)  with  the  information  from  any  one  of 
the  three  points  used  to  evaluate  F  and  M 

p  =  sT(b2  -  a2)/(2a2)  (3.16) 

From  the  above  development,  it  can  be  concluded  that  theoretically 
the  loads  F,  M  and  p  acting  at  the  cross  section  can  be  determined  from 
known  stresses  at  three  points  of  the  external  surface  of  the  cross 
section. 

A  digital  computer  program  called  TIBISAY  was  developed  to  apply 
the  preceding  treatment.   The  program  was  tested  using  a  mathematical 
model  of  a  cantilevered  pipe  as  in  Section  2  of  this  thesis.   Stresses 
generated  at  cross  section  No.  1,  at  points  0°,  120°  and  240°  from  the 
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vertical  downwards  direction  were  used.   The  recovered  loads  were  almost 
identical  with  the  input  ones  producing  the  stresses.   Deviations  are 
due  to  small  round-off  errors  introduced  by  the  computer.   The  program 
and  results  of  the  test  are  presented  in  Appendix  C  of  this  thesis. 
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4.   Conclusions. 

The  results  obtained  in  the  tests  for  applicability  of  the  theories 
developed  in  this  thesis  indicate  that  the  proposed  methods  can  be  used 
for  analysing  stresses  in  pipes. 

There  are,  however,  many  factors  influencing  the  validity  of  the 
results.   Some  of  them  are  related  to  the  amount  of  error  introduced 
when  determining  geometries  and  material  properties.   Others  are  due  to 
proximity  to  concentrated  loads  and  stress  concentrations.   Other  fac- 
tors are  directly  related  to  the  degree  of  redundance  of  strain  gage 
data,  the  number  of  gages  used  and  their  arrangement  around  the  cross 
section,  and  the  number  of  Fourier  coefficients  to  be  used.  Also,  it  may 
be  that  the  least-squares  method  used  does  not  produce  the  best  results; 
perhaps,  by  minimizing  some  other  function  of  error,  better  results  could 
be  obtained. 

First  steps  were  taken  to  investigate  further  some  of  these  factors. 
Appendix  E  of  this  thesis  presents  results  of  a  survey  made  concerning 
the  influence  of  the  number  of  Fourier  coefficients  used  in  the  accuracy 
of  results.   There  are  also  preliminary  results  of  the  influence  of  the 
number  of  gage  elements  and  their  arrangement  around  the  cross  section. 
Due  to  time  limitations,  further  investigations  were  not  conducted. 

The  digital  computer  programs  included  in  this  thesis  can  be  used 
to  investigate  deeper  into  the  field.   Some  refinements  of  the  methods 
can  be  added.   It  is  believed  that  the  following  points  should  not  be 
neglected. : 

a.   Other  criteria  for  minimizing  errors  can  be  investigated. 
A  discussion  of  such  criteria  is  given  in  reference  7  of  this  thesis. 
It  is  likely  that  some  of  the  methods  described  there  may  have  applica- 
tion in  the  pipe  stress  problem  considered  herein. 
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b.  One  may  investigate  the  question  of  the  minimum  number  of 
strain  gage  elements  necessary  to  obtain  acceptable  stress  distributions, 
and  their  orientation  around  the  cross  section.   Similarly,  the  optimum 
number  of  elements  and  preferred  orientation  could  be  investigated. 

c.  While  observing  the  results  for  a  determined  gage  arrange- 
ment, it  was  noticed  that  good  results  were  obtained  using  from  3  to  7 
Fourier  coefficients,  but  more  than  that  produced  bad  results;  using 
again  5  coefficients  and  changing  the  gage  arrangement,  good  results  were 
obtained  for  some  arrangements  and  bad  results  for  others.   Some  of  these 
results  are  presented  in  Appendix  G  of  this  thesis.   Further  study  could 
be  made  of  the  optimal  number  of  Fourier  coefficients  to  be  employed. 

d.  A  method  for  data  rejection  should  be  incorporated  in  the 
theory  and  computer  programs.   There  are  instances  when  it  is  obvious 
that  a  strain  gage  element  has  failed  and  the  reading  of  that  gage  should 
be  neglected.   Usually  bad  gages  are  not  easily  identified  and  bad  data 

is  processed  along  with  good  data.  There  is  the  possibility  of  incorporat- 
ing a  data  rejection  procedure  in  the  method  developed  in  Section  2  of 
this  thesis.   One  suggestion  for  such  a  procedure  is  given  in  the  follow- 
ing. 

The  estimated  Fourier  coefficients  of  the  stress  distributions  are 
evaluated  from  a  minimization  of  the  sum  of  the  squares  of  the  errors, 
i.e.,  the  differences  between  the  gage  element  readings  and  what  they 
should  have  read  theoretically.  We  can  assume  that  these  errors  are 
random  having  a  statistical  normal  distribution  with  mean  and  standard 
deviation  which  could  be  computed.   Establishing  a  level  of  significance 
for  a  deviation  to  be  sufficiently  apart  from  the  mean,  one  can  reject 
the  data  having  greater  deviations  and  recompute  the  stress  distribu- 
tions with  the  remaining  data.   Similar  screenings  can  be  done  with  the 
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newly  computed  stress  distributions  until  the  remaining  data  have  de- 
viations under  the  level  of  rejection.   This  theory  could  be  readily 
used  if  the  redundance  of  data  were  large.   Usually  this  is  not  the  case 
and  one  must  avoid  rejecting  data  which  are  not  perfect  but  which  are 
useful  and  necessary.   If  the  data  is  known  to  have  a  normal  curve  with 
a  very  sharp  peak,  bad  data  is  easily  determined  and  low  significance 
levels  can  be  used.   If  the  bell  is  fairly  flat,  the  majority  of  the 
deviations  have  pronounced  deviations  from  the  mean  and  the  level  of 
rejection  becomes  difficult  to  establish.   The  writer  is  inclined  to 
think  that  there  are  criteria  that  may  satisfy  all  conditions.   Once 
they  are  found,  the  theory  could  be  applied  to  this  problem  and  to 
similar  ones. 
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APPENDIX  A 
DATA  GENERATION  THEORY  AND  RELATED  DIGITAL  COMPUTER  PROGRAMS 
The  data  used  to  test  the  applicability  of  the  theories  developed 
in  this  thesis  was  generated  for  a  mathematical  model  of  a  cantilevered 
pipe  of  concentric  bore  subjected  to  internal  pressure  p  and  general 
force  F  and  general  moment  M  applied  to  the  tip.  Details  of  this  pipe 
are  shown  in  Fig.  A.l.   The  following  material  properties,  pipe  dimensions, 
and  loads  were  specified: 


Fig.  A.l. 
E  =  30,000  ksi 
v  -  0.3 

Specific  weight  -  0.283  lb/in: 
External  diameter  ■  12.0  in 
Internal  diameter  =  10.0  in 


F  =  F  =  F  =  100.0  lb 
x    y    z 

M  =  M  -  M  -  100.0  lb-in 
x    y    z 


p  =  100.0  psi 
Subroutine  CANTIL  presented  herein  in  this  Appendix,  was  used 
to  evaluate  the  tangential,  longitudinal  and  shear  stress  distributions 
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at  the  external  circumference  of  cross  sections  1  and  2  of  this  pipe. 
The  stress  distributions  were  used  in  turn  to  supply  the  data  required 
for  the  tests  specified  in  other  parts  of  this  thesis. 

This  subroutine  is  used  by  either  of  the  (main)  programs  TERESITA 
or  TIBISAY  described  in  Appendices  B  and  C.   The  arrangement  of  the  card 
input  is  indicated  in  Fig.  B.l.  of  this  thesis.   Description  and  specifi- 
cations for  CANTIL  follow  this  page. 
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DESCRIPTION  OF  SUBROUTINE  CANTIL 
Subroutine  CANTIL  evaluates  the  state  of  stress  at  any  point  of 
a  cantilevered  straight  pipe  of  concentric  bore,  with  internal  pressure 
and  loaded  at  the  tip  by  general  forces  and  general  moments.   The  program 
is  capable  of  analysing  from  1  to  360  points  in  each  of  1  to  5  cross 
sections  of  the  pipe.   The  points  have  to  be  at  the  same  radial  distance 
from  the  center  of  the  pipe,  and  separated  from  one  another  by. angular 
intervals  which  are  submultiples  of  360°,  from  1°  to  360°. 

The  data  deck, see  Fig.  B.l, includes  the  following  cards: 
Card  No.  1.  Format  3F20.8 

a.  External  diameter  of  the  pipe,  inches. 

b.  Internal  diameter  of  the  pipe,  inches. 

c.  Pipe  length,  inches. 
Card  No.  2.  Format  3E20.8,  2F10.8 

a.  Young's  modulus  of  elasticity,  psi. 

b.  Shear  modulus  of  elasticity,  psi.   This  can  be 

omitted  in  which  case  it  will  be  computed  by  the 
program  from  E/2(  1  +  v) . 

c.  Coefficient  of  thermal  expansion.   This  can  be  omitted, 
It  is  not  actually  used  in  any  present  application. 

d.  Poisson's  ratio. 

e.  Specific  weight  of  the  material,  pounds  per  cubic 
inch. 

Card  No.  3.   Format  215,  5F10.5 

a.  Number  of  points  to  be  analysed  in  each  cross  section. 

b.  Number  of  cross  sections  to  be  analysed  in  the  pipe. 

c.  The  distances  from  the  tip  of  the  pipe  to  the  cross 
sections  to  be  analysed,  inches. 
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Card  No.  4.   Format  F80.2 

a.   Radial  distance  of  the  points  to  be  analysed,  inches. 
NOTE:  for  applications  to  other  parts  of  this  thesis, 
evaluations  are  made  only  at  the  outside  sur- 
face, but  the  program  is  capable  of  making 
evaluations  at  interior  points. 
Card  No.  5.   Format  F80.2 

a.   Internal  pressure,  psi. 
Card  No.  6.   Format  3F20.8 

a.   Components  of  the  force  applied  at  the  tip  of  the 
pipe,  vertical,  sidewards  and  outwards  from  the 
tip,  so  as  to  form  a  right  handed  set  of  orthogonal 
axes,  lbs. 
Card  No.  7.   Format  3F20.8 

a.   Components  of  the  moment  applied  at  the  tip  of  the  pipe, 
in  the  same  arrangement  as  for  the  forces,  lb-inches. 
Card  No.  8.   Format  F80.2 

a.   This  is  a  sentinel  card  with  one  number  on  it.   A 

negative  number  indicates  that  the  job  is  finished; 
a  zero  or  any  positive  number  indicates  that  the 
program  will  be  rerun  for  another  set  of  radial 
distance  and  loads,  in  which  case,  cards  similar  to 
No.  4  through  8  must  follow  this  one. 
OUTPUTS  (all  arguments  are  outputs) 

RO.  External  radius  of  the  pipe. 
RI.  Internal  radius  of  the  pipe. 
XE.  Young's  modulus  of  elasticity. 
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G.   Shear  modulus  of  elasticity. 
XNU.   Poisson's  ratio. 
SPECW.   Specific  weight  of  the  material. 

D.   A  constant  evaluated  when  analysing  shear  stresses  pro- 
duced by  shearing  forces. 
XL1.   A  one  dimensional  array  containing  the  distances  from 
the  tip  to  each  cross  section  considered. 
NPOINT.   Number  of  points  analysed  in  each  cross  section. 
AREA.   Cross  sectional  area  of  the  pipe. 

ALL1.   A  three  dimensional  array  intended  to  store  the  stresses 
at  any  of  the  points,  in  case  the  points  are  at  the 
external  surface  of  the  pipe;  otherwise,  it  has  no  mean- 
ing. 
Other  subroutines  which  are  needed  and  listings  for  which  are  given  in 
this  Appendix  are: 

TORS:   used  by  CANTIL  to  compute  shear  stresses  produced  by 

torsional  moments. 
BEND:   used  by  CANTIL  to  compute  bending  stresses. 
AXIAL:   used  by  CANTIL  to  compute  stresses  due  to  axial  forces. 
SHEAR:   used  by  CANTIL  to  compute  shear  stresses  produced  by 

shear  forces. 
CROSS:   computes  the  cross  product  of  two  vectors. 

DOT:   computes  the  dot  product  of  two  vectors. 
TIMES:   computes  the  product  of  a  scalar  and  a  vector. 

PLUS:   computes  the  sum  of  two  vectors. 
MINUS:   computes  the  difference  between  two  vectors. 
TRANSF:   computes  the  components  of  a  vector  in  another  angularly 
displaced  set  of  axes. 
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METHOD. 

The  vectorial  approach  discussed  in  reference  5  of  this  thesis  is 
used  to  compute  stresses  produced  at  points  of  the  cross  section  by 
internal  pressure  and  forces  and  moments  acting  at  the  centroid  of  that 
cross  section. 


-A 


(a) 


-f; 


1 

/ 


(b) 


-M, 


Fig.  A-2. 

Those  forces  and  moments,  in  turn,  are  computed  by  the  classical  vector!' 

al  analysis  of  the  statics  of  a  beam.   If  the  force  F  and  the  moment  M 

are  applied  at  one  end  of  a  beam,  Fig.  A. 2. a.,  the  force  F  and  moment 

n 

*M  acting  at  the  centroid  of  the  cross  section  nn,  Fig.  A.2.C.,  are 
n 

found  by  considering  the  equilibrium  of  the  segment  of  beam  from  the 
tip  to  the  cross  section.   Summing  forces  and  taking  moments  about  the 
centroid  of  the  cross  section 
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F  +  W  +  Fn  =  0  (A.l) 

M   +  M  +  LxF  +  %LxW  =  0  (A.  2) 

n 

where  L  is  the  vector  distance  from  the  centroid  of  the  cross  section 

to  the  point  of  application  of  the  loads  at  the  tip,  and  W  is  the  force 

due  to  the  weight  of  the  segment  of  beam. 

Solving  for  F  and  M 
°      n      n 

F   -  -  (F  +  W)  (A. 3) 

M  -   -  (M  +  LxF  +  %LxW)  (A. 4) 

The  forces  and  moments  acting  at  the  cross  section  under  considera- 
tion on  the  remaining  portion  of  the  beam,  are  the  negatives  of  F  and 

M  found  above,  Fig.  A.2.b. 
n 

PRINT -OUTS 

The  print-outs  are  several  pages  as  follows 
1st  page. 

a.  Pipe  dimensions 

b.  Computed  cross  sectional  area,  moment  of  inertia 
and  polar  moment  of  inertia  of  the  cross  section. 

c.  Material  properties, 

d.  Number  of  points  analysed  in  each  cross  section, 
number  of  cross  sections  analysed  and  their  distances 
from  the  tip. 

2nd  page. 

a.  Radial  distance  of  the  points  analysed  . 

b.  Loads  applied  to  the  tip. 

c.  Distance  from  the  tip  to  the  first  cross  section 
analysed. 
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d.   Computed  loads  acting  at  that  cross  section 
Following  pages.   For  each  analysed  point  the  print-out  contains: 

a.  Angular  position  of  the  point  in  the  cross  section 
(degrees) 

b.  Forces  and  moments  acting  at  that  cross  section  in 
radial,  tangential  and  axial  components. 

c.  Stresses  produced  by  individual  loads. 

SHTZT  is  the  shear  stress  in  the  ZT  direction 
produced  by  the  axial  component  of  the  moment 
SBZZ  axial  (ZZ  direction)  stress  produced  by 
the  tangential  component  of  the  moment. 
SAZZ  axial  stress  produced  by  the  axial  compon- 
ent of  the  force. 

SHFZR  and  SHFZT  are  stresses  produced  by  the 
radial  and  tangential  components  of  the  force. 
SPRR  and  SPTT  are  radial  and  tangential  stresses 
produced  by  the  internal  pressure. 

d.  If  the  analysed  points  are  in  the  external  surface 
of  the  pipe,  there  will  be  additional  print-outs  of 
the  total  axial,  tangential  and  shear  stresses  at 
the  point. 

Analysis  of  other  cross  sections  for  the  initial  and/or  other  radial 
distance  and  tip  loads  and  pressure  produce  similar  outputs  as  described 
from  page  2  on.   The  complete  program  and  sample  outputs  of  the  analysis 
made  for  points  at  the  external  surface  and  %  in.  under  that  surface  of 
cross  sections  1  and  2  of  the  pipe  shown  in  Fig.  A.I.,  follow  this  page. 
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EXTERNAL 
INTERNAL 
LENGTH 


DIAMETER 
DIAMETER 


12.00  INCHES 
10.00  INCHES 
90.00  INCHES 


CRUSSFCTIONAL  AREA  = 
MOMENT  JF  INERTIA 
POLAR  MOMENT  OF  I.  = 


34.5575  IN2 

527.0020  IN4 

1054.0039  IN4 


MODULUS  OF  ELASTICITY      0.30E  08  PSI 

SHEAR  MODULUS  OF  ELASTICITY      P.12F  08  PSI 

COEFFICIENT  OF  THERMAL  EXPANSION  NOT  GIVEN 

POISSON  ^ATIO     0.300 

SPECIFIC    WFIGHT  0.2R3C    LB/IN3 


NUMBER  OF  QUINTS  TO  ANALYSF 
NUMBER  OF  SECTIONS  TO  ANALYZE 
DISTANCES  FRCM  THE  TIP  TO  THE 


30. OOP   60. 0^0 


0.0 


12 
2 

CROSS 

CO 


SECTIONS,  INCHES 
CO 


53 


RADIAL  DISTANCE  OF  POINTS  = 


6.00000  INCHES 


LOADS  APPLIED  TO  THE  TIP 

INTFRNAL  PRESSURE   = 
APPLIED  FORCE 

APPLIED  MOMENT 


100.00  PSI 

100.00 

100.00 


100.00 
100. 00 


100.00  LB 
100.00  LB-IN 


DISTANCE  FROM  THE  TIP  TO  THE  CROSS  SECTION 


30. CO  INCHFS 


LOADS  ACTING  ON  THE  CROSS  SECTION 
INTERNAL  PRESSURE   =       100.00  PSI 
APPLIED  FORCE       =       393.39       100.00 
APPLIED  MOMENT      =     -2900.00      7500.89 


100.00  LB 
100.00  LB-IN 
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THTT\     = 
1  .      \SI 
1  k  V-   >F.j 

SHTZ1 
S    i  11 

S  '.// 
bHI  ?f 
SHf  /T 

S    r  .  j 

S    L  T  T 


stf  f-  ^s    ; 


'      •     hf     CE 


DFGf rrs 
=    (  .  39339E    '   3 
=-0.29C( OE    04 


.  >t>926F 
r .S5  399F 

.23 3 3 7E 
O.P 
) .  5  3  8  M  E 


C-2 

i   1 

01 


.  * 545  5!     C  2 

T  •  T  PI  '•  h  C  T  I  UN 
ALl  1  (  1  t  1,1) 
? . 4  5  45  5  E      )  3 


o.  iooooe 

0.75009E 


03 


7'Z    DIRECTION 
ALLK2,     1,1) 

-0.825C5E    ">  2 


looo^e 

LOOOOE 


0  3 

03 


LB 

LD-IN 


Z'T    DIRECTION 
ALLK3,     1,1) 

0.59553E    01 


1  HCT  a    =  30.00         DEGF rES 

I        ■   ;l        •  i    :    f-ORCE       =    C  .  39069E 

T-1  IN  >l    UM(    •      ICMENT    =    C . 12  390E 


03-0.11009F    ^3 
04    0. 79460  E    "A 


O.IOOO^F    03    LB 
C.IjOOOE    0  3    LB- IN 


SHTZT 

S  :/  ' 
<-'  A  /  7 
SHF  7  - 
SHf  ZT 
S  ' 
S     n  TT 


r>  h  9  2  6  F 
2893 7F 
5929 7F 
4545 5E 


r  r 

0  2 

ri 

c  1 


ST'  f  S  S ! 


T«T    DIRFCTION 
\  L  L  1  (  1  ,     2,1) 

->,  ,AS<+SSr-     33 


Z'Z    DIRECTION 
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THtTA  =   120.00    OEGPEES 

TRANSFORMED  FORCE   =-C.HCC9E  03-0.39069F  03  C.10000E  03  LB 

TRANSFORMED  MOMENT  =  r.794bOE  04-0.12390E  ^4  0.10000F  03  LR-IN 
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RADIAL  DISTANCE  OF  POINTS  =      6.00000  INCHFS 


LOADS  APPLIED  TO  THE  TIP 

INTERNAL  PRESSURE   =       100.00  PSI 

APPLIED  FORCE       =       LOO. 00       100.00       100.00  LB 

APPLIED  MOMENT      =       100.00       100.00       100.00  LB-IN 


DISTANCE  FROM  THE  TIP  TO  THE  CROSS  SECTION      60.00  INCHES 

LOADS  ACTING  ON  THE  CROSS  SECTION 

INTERNAL  PRESSURE   =       100.00  PSI 

APPLIED  FORCE       =       686.79       100.00       100.00  LB 

APPLIED  MOMENT      =     -5900.00     23703.59       100.00  LB-IN 
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THETA  =     O.n     DFGREBS 

TRANSFORMED  FORCE   =  C68679E  03  O.IOOOOE  03  C.10000E  03  LB 

TRANSFORMED  MOMENT  =-C.59CC0E  ^4  0.23704E  05  0.10000E  03  LB-IN 
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RADIAL  DISTANCE  OF  POINTS  =      5.50000  INCHES 


LOADS  APPLIED  TO  THE  TIP 


INTFRNAL  PRESSURE   =       100.00  PSI 

APPLIED  FORCE       =       100.00       100.00       100.00  LB 

APPLIED  MOMENT      =       100.00       100.00       100.00  LB-IN 


DISTANCE  FROM  THE  TIP  TO  THE  CROSS  SECTION      30.00  INCHES 


LOADS  ACTING  ON  THE  CROSS  SECTION 

INTERNAL  PRESSURE   =       100.00  PS! 

APPLIED  FORCE       =       393.39       100.00       100.00  LB 

APPLIED  MOMENL =     -2900.00      7500.89       100.00  LB-IN 
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RADIAL  DISTANCE  OF  POINTS  =      5.50000  INCHES 


LOADS  APPLIED  TO  THE  TIP 


INTERNAL  PRESSURE   =       100.00  PSI 

APPLIED  FORCE       =       100.00       100.00       100.00  LB 

APPLIED  MOMENT      *       100.00       100.00       100.00  LB-IN 


DISTANCE  FROM  THE  TIP  TO  THE  CROSS  SECTION      60.00  INCHES 

LOADS  ACTING  ON  THE  CROSS  SECTION 

INTERNAL  PRESSURE   =       100.00  PSI 

APPLIED  FORCE       -       686.79       100.00       100.00  LB 

APPLIED  MOMENT      =     -5900.00     23703.59       100.00  LB-IN 
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TH^TA    =       120.00  DEGREES 

TRANSFORMED    f-OPCE       =-0.25679F 
TRANSFORMED    Mi'HENT    =    0.23478E 
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APPENDIX  B 

DIGITAL  COMPUTER  PROGRAM  FOR  SECTION  2  AND  NUMERICAL  RESULTS 

The  programs  presented  herein  provide  a  method  for  applying  the 
theory  developed  in  Section  2  of  this  thesis.   Actual  situations  can  be 
handled  by  subroutines  ZULIA,  GUAYRA,  CUMANA  and  DSIMQ  (DSIMQ  is  present- 
ed in  Appendix  D) .   Simulated  situations  require  in  addition  subroutines 
CANTIL  presented  in  Appendix  A,  and  TIUNA  and  BLANCA  presented  in  this 
Appendix.   The  last  of  the  subroutines  mentioned  is  an  adaptation  of  a 
subroutine  taken  from  reference  7  of  this  thesis. 

ZULIA  is  used  to  specify  the  required  information,  such  as  material 
properties,  number  of  Fourier  coefficients  to  be  computed  and  descrip- 
tion of  the  arrangement  of  strain  gage  elements  and  their  transverse 
sensitivity  factor.   Enough  self  explanatory  comment  statements  have  been 
included  which  specify  the  introduction  of  the  information.   Print-out 
statements  present  the  information  received  for  processing.   Once  the  in- 
formation is  properly  arranged,  subroutine  GUAYRA  is  used  to  make  the  re- 
quired computations. 

The  sequence  in  GUAYRA  follows  closely  the  theoretical  development 
of  the  method.   It  has  sufficient  comment  statements  separating  blocks 
of  computer  processing  instructions  as  to  follow  the  theoretical  pro- 
cedure.  Print-out  statements  present  the  computed  Fourier  coefficients. 

CUMANA  is  used  by  GUAYRA  for  auxiliary  computations. 

Use  of  the  subroutines  is  made  by  a  very  brief  main  program  TERESITA. 
This  program  merely  dimensions  the  various  subroutines  and  turns  control 
over  to  subroutine  ZULIA.   The  card  arrangement  of  the  source  deck  is  as 
shown  in  Fig.  B. 1. 

Some  sample  results  of  the  application  of  the  program  as  specified 
at  the  end  of  Section  2,  follow  the  card  listings  of  the  programs. 
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THETA  = 
STRFSSES 
ACTUAL 
COMPUTED 


DFGREFS 

TANG. 
0.45455E 
0.45455E 


POINT  NO.    1 

LONG. 
C3       0.17552E  01 
C3       0.17554E  01 


SHEAR 
0.59553E    01 
0.59553E    01 


THETA  = 

STRESSES 

ACTUAL 

COMPUTED 


10  DEGREES 

TANG. 
0.45455E 
0.45455E 


PCINT  NO. 


C3 

C3 


11 

LONG. 
0.19702E 
0. 19704E 


01 
01 


SHEAR 
0.49382E    01 
0.49382E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


20  DEGREES 

TANG. 
0.45455E 
0.45455E 


PCINT  NO. 


C3 
C3 


21 

LONG. 

0.22133E  01 

0.22135E  01 


SHEAR 
0.37884E    01 
0.37884E    01 


THETA  = 

STRESSES 

ACTUAL 

COMPUTED 


30  DFGRFFS 

TANG. 
0.454*5E 
0.45454E 


POINT    NO. 


C3 
C3 


31 
LONG. 
0.24770E    01 
0.24773F    01 


SHEAR 
0.25407E 
0.25407E 


01 
01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


40  DEGREES 

TANG. 
0.4*455E 
0.45454E 


PCINT  NO. 


C3 
C3 


41 
LONG. 
0.27534E  01 
0.27538E  01 


SHEAR 
123 3  IE 
12331E 


01 
01 


THETA  = 

STRESSES 

ACTUAL 

COMPUTED 


50  DEGREES 

TANG. 
0.45455E 
0.45454E 


PCINT  NO. 


C3 

C3 


51 
LONG. 
30341E  01 
30345F  01 


SHEAR 
0.946 09 E-01 
0.94599E-01 


THETA  = 

STRESSES 

ACTUAL 

COMPUTED 


60  DEGREES 

TANG. 
0.45455E 
0.45455E 


PCINT  NO. 


C3 
C3 


61 

LONG. 
33104F  01 
33110E  01 


SHEAR 
0.14022E    01 
0.14.022E    01 


THETA  - 
STRESSES 
ACTUAL 
COMPUTED 


70  DEGRFFS 

TANG. 
0.45455E 
0.45455E 


POINT    NO. 


C3 
C3 


71 
LONG. 
0.35742E    01 
0.35748E    01 


SHEAR 
0.26498F 

■0.26498E 


01 
01 


THETA  = 
STRESSES 
ACTUAL 
CCMPUTED 


80  DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 
C3 


81 
LONG. 
0.38172E 
0.38178E 


01 
01 


SHEAR 
•0.37997E    01 
0.37997E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


90  CEGREFS 

TANG. 
0.45455E 
0.45455E 


PCINT    NO. 


C3 
C3 


91 

LONG. 
0.40322E 
0.40328E 


01 
01 


SHEAR 

■0.48168E 

0.48168E 


01 
01 


THETA  =  ICO 
STRESSFS 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45455E  C3 


POINT  NO.  101 

LONG. 
0.42126E  01 
0.42132E  01 


SHEAR 
0.56702E 
0.56702E 


01 

01 


THETA  =  110 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45455E  C3 


POINT  NO.  Ill 

LONG. 
0.43530E  01 
0.43535E  01 


SHEAR 

0.6  33  41  E 

•0.63341F 


01 
01 
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THETA  =  120  DEGREES  POINT  NO.  121 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.44490E  01      -0.67882E  01 

COMPUTED  0.45455E  C3  0.44495E  01      -0.67882E  01 

THETA  =  130  DEGREES  POINT  NO.  131 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.44977E  01      -0.70188E  01 

COMPUTED  0.45455E  C3  0.44982E  01      -0.70188E  01 

THETA  =  140  DEGREES  POINT  NO.  141 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.44977E  01      -0.7C188E  01 

COMPUTED  0.45455E  C3  0.44983E  01      -0.70188E  01 

THETA  =  150  CEGREES  POINT  NO.  151 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.44490E  01      -0.67882E  01 

COMPUTED  0.45455E  C3  0.44496E  01      -0.67882E  01 

THETA  =  160  DEGREES  PCINT  NO.  161 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.43530E  01      -0.63341E  01 

COMPUTED  0.45455E  C3  0.43537E  01      -0.63341E  01 

THETA  =  170  DEGREES  POINT  NO.  171 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.42126E  01      -0.56703E  01 

COMPUTED  0.45455E  C3  0.42133E  01      -0.56702E  01 

THETA  =  180  DEGREES  POINT  NO.  181 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.40322E  01      -0.48168E  01 

COMPUTED  0.45455E  C3  0.40328E  01      -0.48168E  01 

THETA  =  190  DEGREES  POINT  NO.  191 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.38172E  01      -0.37997E  01 

COMPUTED  0.45455E  C3  0.38178E  01      -0.37997E  01 

THETA  =  200  CEGREES  POINT  NO.  201 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.35742E  01      -0.26498E  01 

COMPUTED  0.45455E  C3  0.35747E  01      -0.26498E  01 

THETA  =  210  DEGREES  POINT  NO.  211 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  C.45455E  C3  0.33105E  01      -0.14022E  01 

COMPUTED  0.45455E  C3  0.33110E  01      -0.14022E  01 

THETA  =  220  CEGREES  POINT  NO.  221 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.30341E  01      -0.94623E-01 

COMPUTED  0.45455E  C3  0.30347E  01      -0.94618E-01 

THETA  =  230  CEGREES  POINT  NO.  231 

STRESSES  TANG.             LONG.             SHEAR 

ACTUAL  0.45455E  C3  0.27534E  01       0.12331E  01 

CCMPUTED  0.45455E  C3  0.27541E  01       0.12331E  01 
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THETA  =  240 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO.  241 

LONG. 
C3       0.24770E  01 
C3       0.24777E  01 


SHEAR 
0.25407E 
0.25407F 


01 
01 


THETA  =  250 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E 
0.45455F 


POINT    NO 


C3 
C3 


251 

LONG. 
22133E    01 
22139F    01 


SHEAR 
0.37883E    01 
0.37R81E    01 


THETA  =  260 

STRESSES 

ACTUAL 

COMPUTED 


CEGREES 

TANG. 
0.45455E 
0.45455E 


PCINT    NO. 


C3 
C3 


261 

LONG. 
19702F    01 
19708E    01 


SHEAR 
0.49382E 
0.49382E 


01 
01 


THETA  =  270 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 

0.45455E 


POINT  NO.  271 

LONG. 
C3       0.17552E  01 
C3       0.17558E  01 


SHFAR 
0.59553F 
0.5Q553E 


01 
01 


THETA    =    280 

STRESSES 

ACTUAL 

COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 
C3 


281 

LONG. 
15748E  01 
15754F  01 


SHEAR 
0.68088F 
0.68087F 


^1 
01 


THETA  =  290 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E  C3 
0.45455E  C3 


POINT    NO.     291 

LONG. 
0.14345E    01 
0.14351E    01 


SHEAR 
0.74726E 
0.74726E 


01 
01 


THETA  =  3C0 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45455E  C3 


POINT  NO.  301 

LONG. 
0.13385E  01 
0.13392E  01 


SHFAR 
0.79267E 
0.79267E 


01 
01 


THETA  =  310 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT    NO.     311 

LONG. 
C3  0.12898E    01 

C3  0.12905E    01 


SHEAR 
0.81573F    01 
0.81573E    01 


THETA  =  320 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E  C3 
0.45455E  C3 


POINT  NO.  321 

LONG. 
0.12897E  01 
0.12904E  01 


SHEAR 
0.81573F 
0.815 73F 


01 

01 


THETA  =  330 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45455E  C3 


POINT  NO.  331 

LONG. 
0.13385E  01 
0.13390F  01 


SHEAR 
0.79267E    01 
0.79267E    01 


THETA  =  340 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E 
0.45455F 


POINT  NO. 


C3 
C3 


341 

LONG. 
14345E  01 
14348E  01 


SHEAR 
0.74726E    01 
0.74726E    01 


THETA  =  350 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E  C3 
0.4545^E  C3 


POINT  NO.  351 

LONG. 
0.15748F  01 
0.15751E  01 


SHFAR 
0.680 8RE 
0.68088E 


01 

01 
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THFTA  = 
STRESSES 
ACTUAL 
COMPUTED 


0  DEGREES 

TANG. 
0.4S4^6F 
0.45455F 


POINT  NU 


0  3 
03 


1 
LONG. 
17552E  01 
L7556E  01 


SHEAR 
0.59553E  01 
0.59553E  01 


THt TA  =   U 
STPEbSEo 
ACTUAL 
COMPUTED 


CHOREES 

TANG. 
0.4S455E 

0.45455E 


POINT    NO 


03 

03 


11 

LONG. 
1^>702E    CI 
1Q706E    01 


SHEAR 
<^.49382E 
0.49382E 


01 
01 


THtrTA     = 
STRESSES 
ACTUAL 
COMPUTED 


20  DEGREES 

TANG. 
0»45455E 

0  .  4  5  4  5  5  E 


POINT  NO 


03 

03 


21 

LONG. 
22133E  01 
22137F  01 


SHEAR 
C.37884E 
0.37884F 


01 
01 


THFTA  = 
STRESSES 
ACTUAL 
COMPUTE  I) 


3C  OF 


GREES 

TANG. 
0.4S455L 
0.45455F 


POINT  NO 


C3 
03 


31 

LONG. 

0.24770E  01 

0.24774E  01 


SHEAR 
0.2^407E 
C.25407E 


01 

01 


THFTA  = 
STRESSES 
ACTUAL 
COMPUTED 


4  0    DEGREES 

TANG. 
0.454S5E 
0.45455E 


POINT  NO 


03 
0  3 


41 

LONG. 
27534F  01 
27538E  01 


SHEAR 
0.12331E  01 
0.1 233  IE  01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.4545^ 
0.4545  5E 


POINT  NO 


03 

03 


51 

LONG. 
30341E  01 

30345E  01 


SHEAR 

•0.9460^-01 

0.94604E-01 


THFTA  = 
STRESSES 
ACTUAL 
COMPUTED 


60  DFGREES 

TANG. 
0.45455E 
0.45455F 


POINT  NO. 


03 
03 


61 

LONG. 
33104E  01 

33109E  01 


SHEAR 

0.14C22E    01 
0.14022E    01 


THFTa  = 
STRESSES 
ACTUAL 
COMPUTED 


70  CFoPFES 

TANG. 
C.4^45*>E 

0.45455F 


POINT  NU 


03 

03 


71 
LONG. 
35742E  01 
35746E  01 


SHEAR 
•0.26498E 
0.26498E 


01 
01 


THFTA  = 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.4545'JE 

0.45455F 


POINT  NU 


U3 
03 


81 

LONG. 

0.38172E  01 

0.38177E  01 


SHEAR 
•G.37997E 
■0.37997E 


01 
01 


THFTa  = 
STRESSES 
ACTUAL 
COMPUTE,) 


9C    DEGREES 

TANG. 
0.45455E 
G.45455E 


POINT  NO 


03 
03 


91 
LONG. 
40322E  01 
40327E  01 


SHEAR 
•D.48168E 
•0.48168E 


01 

01 


THETA  =  IOC 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
G.45455E 
0.45455F 


POINT  NO 


03 

03 


101 

LONG. 
0.42126E  01 
0.42132E  CI 


SHEAR 

■0.56702E 

0.567C2E 


01 
01 


THETA  =  110 
STRESSES 
ACTUAL 
COMPUTED 


DEGREE.S 

TANG. 
0.45455F 
0.45455F 


POINT  NO 


03 
03 


111 

LONG. 
43530E  01 
43535E  01 


SHEAR 
0.63341E  01 
0.63341F  01 
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THETA  =  120 
STRESSES 
ACTUAL 
COMPUTET 

LEGREES 

TANG. 
U.45455E 
0.45455F 

POINT 

03 
03 

NO.  121 

LONG. 
0.44490E 
0.44495F 

01 
01 

SHEAR 
-0.67882E 
-P.67R82E 

01 
01 

THETA  =  130 

STRESSES 

ACTUAL 

COMPUTE! 

DEGREES 

TANG. 
0.45455E 
0.45455E 

POINT 

03 
03 

NO.  131 

LONG. 
0.44977E 
0.449R3E 

01 
01 

SHEAR 
-0.70188E 
-0.70188F 

01 
01 

THETA  =  140 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455E 

POINT 

03 
03 

NO.  141 

LONG. 
0.449776 
0.44983E 

01 
01 

SHEAR 
-0.70188E 
-0.70188E 

01 
Oi 

THETA  =  150 

STRESSES 

ACTUAL 

COMPUTE D 

DEGREES 

TANG. 
0.45455F 
0.45455F 

POINT 

03 
03 

NO.  151 

LONG. 
0.44490E 
0.44496E 

01 
CI 

SHEAR 
-0.67882E 
-0.67882E 

01 
01 

THETA  =  160 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455E 

POINT 

03 
03 

NO.  161 

LONG. 
0.43530E 
0.43536E 

01 
CI 

SHEAR 
-0.63341E 
-0.63341E 

01 
01 

THFTA  =  170 
STRESSES 
ACTUAL 
COMPUTED 

CFGREES 

TANG. 
0.<t5455F 
0.45455F 

POINT 

03 
03 

NO.  171 

LONG. 
0.42126E 
0.42133E 

01 

CI 

SHEAR 
-0.56703E 
-0.56702E 

01 
01 

THETA  =  180 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455E 

POINT 

03 
03 

NO.  181 

LONG. 

0.40322E 

.  0.40329E 

01 
CI 

SHEAR 
-0.48168E 
-0.48168E 

01 
01 

THETA  =  190 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455E 

POINT 

03 

03 

NO.  191 

LONG. 
0.38172E 
0.38179E 

01 
01 

SHEAR 
-0.37997E 
-0.37997E 

01 
01 

THETA  =  200 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455F 

POINT 

03 

03 

NO.  201 

LONG. 
0.35742E 
0.35749E 

01 

01 

SHEAR 
-0.26498F 
-0.26498E 

01 
01 

THETA  =  210 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45456E 
0.454«S5E 

POINT 

03 
03 

NO.  211 

LONG. 
0.33105E 
0.33111E 

01 
01 

SHEAR 
-0.14022E 
-0.14022E 

01 

01 

THETA  =  220 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455E 

PL)  I  NT 

03 
03 

NO.  221 

LONG. 
0.3C341E 
0.30347E 

01 

01 

SHEAR 
-0.94623E- 
-0.94625E- 

oi 
01 

THETA  =  230 
STRESSES 
ACTUAL 
COMPUTED 

DEGREES 

TANG. 
0.45455E 
0.45455E 

POINT 

03 
03 

NO.  231 

LONG. 
0.27534E 
0.27541F 

01 
01 

SHEAR 
0.12331E 
0.L2331E 

01 
01 
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THt-TA  =  240 
STRESSES 
ACTUAL 
COMPUTEJ 


DEGREES 

TANG. 
0.45455F 
0.45455F 


POINT  NO 


03 
03 


0 


241 

LONG. 
24770E  01 
24777F  01 


SHEAR 
0.25407E 
0.25407F 


01 
01 


1HET\  =  2  50 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 
0  . 45455L 


POINT    NO 


03 
02> 


251 

LONG. 
22133E  01 
22139E  CI 


SHEAR 
0.37883F 
0.37883E 


01 
01 


IHETA     =     ; 
STRESSES 

ACTUAL 
COMPUTFD 


160  DFGREES 

TANG. 
0.45455E 

0.45455F 


POINT  NO 


03 

03 


261 

LONG. 
19702F  01 
19708E  01 


SHEAR 
0.49382F 
0.49382F 


01 
01 


THETA  =  2  70 
STRESSES 
ACTUAL 
COMPUTED 


DFGREES 

TANG. 
0.45455E 
0.454551 


POINT  NO 


03 
03 


271 

LONG. 
17552E  01 
17558F  01 


SHEAR 
0.59553E  01 
0.59553E  01 


THt-T«  =  <>*> 
STPFSSFS 
ACTUAL 
COMPUTE  ) 


DFGREES 

TANG. 
0.45455E 
0.454b5F 


POINT  NO 


03 
03 


281 

LONG. 
15748F  01 
15754E  01 


SHEAR 
0.68088E  01 
0.68087E  *l 


THbTA  =  290 
SU  ESSES 
ACTUAL 
COMPUTED 


CFGRFES 

TANG. 
0. 454S5F 
0.4  54  5  5E 


POINT  NO.  291 


03 
03 


LONG. 
14345F  01 
1435CE  01 


SHEAR 
0.74726E 
0.74726F 


01 

01 


7HRTA  =     300 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455F 
0.45455F 


POINT  NO. 


03 
03 


301 

LONG. 
13385E  01 
13390E  01 


SHEAR 
0.79267E  01 
0.79267E  01 


THLTA  =  3  10 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455F 
0.45455E 


POINT  NO 


03 
03 


311 

LONG 
.28981 
12903E  01 


SHEAR 
0.81573E  01 
0.81573E  01 


THETA  =  320 
STRESSES 
ACTUAL 
COMPUTED 


UFGRFFS 

TANG. 
0.454S5E 
0.45455F 


POINT  NO. 


03 
03 


321 

LONG. 
L2897E  CI 
12902F  01 


SHEAR 
0.81573E  01 
0.81573E  01 


THETA  =  330 
STRESSFS 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.4(3455F 
0.454r>5F 


POINT  NO 


03 
03 


331 

LONG. 
13  38  5E  01 
13390E  01 


SHEAR 
0.79267F 
0.79267F 


01 


THFTA  =  340 
STRESSES 
ACTUAL 
COMPUTED 


iEGREES 

TANG. 
0.*5455E 
0.45455F 


POINT  NO 


03 
03 


341 

LONG. 
14345E  01 
14349E  01 


SHEAR 
0.74726F  01 
0.74726F  01 


THETA  =  3  50 
STRESSES 
ACTUAL 
COMPUTE  ") 


DEGREES 

TANG. 
0.45455F 
0.45455F 


POINT  NO 


03 
03 


351 

LONG. 
15748E  01 
15752E  01 


SHEAR 
0.68088E  01 
C.68088E  01 
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THtTA  = 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455F 
0.45454E 


POINT  NO. 


03 
03 


1 

LONG. 

0.17552E  01 

0.17556E  01 


SHEAR 
0.59553E  01 
0.59553E  01 


THFTA  = 
STRESSES 

ACTUAL 
COMPUTE  ) 


10  DEGREES 

TANG. 

.45455E 

C  .45454E 


POINT  NO.   11 

LONG, 
03       C.19702E  CI 
03       0.19706E  01 


SHEAR 
0.49382E  01 
0.49  382E  01 


THFTA  =  2C  CEGRFLS  POINT  NO.   21 

STRFSStS  TANG.  LONG. 

ACTUAL  1.45455E  03       0.22133E  01 

COMPUTED  r.454S4F  03       0.22136E  CI 


SHEAR 
0.37884E 
0.37K84E 


01 

01 


THFTA  = 
STRESSES 

ACTUAL 
COMPUTED 


30  DEGREES 

TANG. 
G.45455E 
0.45454F 


POINT  NO.   31 

LONG. 
03       0.24770E  01 
03       0.24774E  01 


SHEAR 
0.25407E 
0.25407E 


01 
01 


THtTA  = 
STRESStS 
ACTUAL 
COMPUTED 


FGREES 

TANG. 
0.45455F 
D.45454E 


PUINT  NO.   41 

LONG. 
03       0.27534E  01 


03 


0.27538E  01 


SHFAR 
0.12331E  01 
0.12331E  01 


THETA  = 
STFESSES 
ACTUAL 
COMPUTED 


50  DEGREES 

TANG. 
0.45455E 
C  .45454E 


POINT  NO.   51 

LONG. 
03       0.30341E  01 
03       0.30344E  01 


SHEAR 

0.94609F-01 

•0.94596E-01 


THFTA  = 
STRESSES 
ACTUAL 
COMPUTED 


60  DEGREES 

TANG. 
0.4 54 5 5F 
C .45454E 


POINT  NO 


03 
03 


61 
LONG. 
33104E  01 
33108E  01 


SHEAR 

•0.14022E    01 

0.14022E    01 


THFTA  = 
STRFSSES 
ACTUAL 
COMPUTED 


7C     DEGREES 

TANG. 
0.45455F 
D.45454F 


POINT    NO 


03 
03 


71 

LONG. 
35742E    CI 
35746E    01 


SHEAR 

■0.26498E    ^1 

C.26498E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTE ) 


80  DEGREES 

TANG. 
0.45455F 
0.45454E 


POINT    NO. 


03 

03 


81 

LONG. 
0.38172E    01 
0.38177E    01 


SHEAR 

0.37997E    01 

■0.37997E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTE  ) 


90  DEGREES 

TANG. 
0.45455E 
0.45454F 


POINT    NO. 


03 
03 


91 
LONG. 
40322E    01 
4C327E    01 


SHFAR 
■0.48168E 
0.48168E 


01 
01 


THETA  =  IOC- 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
U.45454E  03 


POINT  NO.  101 

LONG. 
0.42126E  01 
0.42131F  01 


SHEAR 
0.56702E  01 
0.56703E  01 


THETA  =  11 
STRESSES 
ACTUAL 
COMPUTE  1 


DEGREES 

TANG. 
0.45455F 
0.45455F 


POINT  NO.  Ill 

LONG. 
03       0.43  5  30E  0  1 
03       0.43535E  01 


SHEAR 
0.63341E  01 
0.63341E  01 


95 


mm  roomer 
oooooo 

LLLU  LULL  LULU 

oooooo 
oooooo 
oooooo 

OOlf^r-m 
mtr,«-«M(\im 

•  ••••• 

oooooo 


f  mmmmm 

oooooo 

LULL'LUU-LUU 

oooooo 
oooooo 
oooooo 

<\J<\I(\IC\I(\J(\! 


•  ••••• 

oooooo 


*C  *£■  *i0  nC  ^C1  *£ 

oooooo 
I  I  I  I  I  I 

U-UJUJILUJUJ 

co«$-un«©om 

O0D.-<\©Sf-m 

mmmcr^"-' 

m>cr»->$,i«-'t-i 

•  ••••• 

oooooo 


MMmmmm 

MMfMfMfVjrvj 

*0<C  *o»o^«o 

oooooo 

oooooo 

oooooo 
1  1  1  1  1  1 

LULU  LULU  LULU 

LLLU  LULU  LULU 

LULULU  LULULU 

oooooo 

oooooo 

>t>*mr-oco 

oooooo 

oooooo 

OMO-<0CO 

oooooo 

oooooo 

C^MIXN^-O^OO 

ooLo«-»^m 

oooooo 

mr>r\i©^t«-< 

mc>'-<M(\jm 

s£)  sQ  sC  *G  >£)  vC 

lTiMM>frsCr» 

•  ••••• 

•    ••••• 

•    ••••• 

LL.1 

oooooo 

oooooo 

oooooo 

CO 

3E 

i 

<M<\immmm 

<*>*>*>**■>* 

JC 

O 

0C 

oooooo 

oooooo 
1  1  1  1  1  1 

LULU  LLLU  LULU 

LULU  LULU  LULU 

Q 

W 

o 

OOOOOO 

r^vOmo^r*- 

C 

OOOOOO 

ON«00(T(\1 

< 

oooooo 

0*  r-»aooo©m 

UJ 

OOiPHsen 

•4"  ITi  >*•  <*  l/MT\ 

ac 

mc^— <rvif\im 

oooooo 

g 

i— <  ^-  .—  — <  -H  --4 

•  ••••• 

•    ••••• 

•    ••••• 

n> 

z 

oooooo 

oooooo 

oooooo 

o 

UJ 

o 

K 

vo 

o 

csjmmmm 

— 

mm mm mm 

• 

^  *c^o>&  *b*o 

N- 

H 

UJ 

ooooo 

X 

oooooo 

UJ 

oooooo 

co 

CO 

CO 

1  1  1  1  1  1 

a 

CO 

LU  LULU  LULU 

-x 

LULULU  LULU  LL 

»— 1 

UJ  LULU  LULULU 

UJ 

CO 

ooooo 

o 

oooooo 

z 

v©v0if\vOCJy(\l 

00 

o>o 

CO 

ooooo 

OC 

oooooo 

oor^cr^i—io 

o 

<im 

C 

ooooo 

oooooo 

2 

ir»rsi<— «\J>t"M 

o 

or 

OMCO^O 

o 

cgr\ir\!r\jr\jr\i 

o 

o^om^m** 

i-iUJ 

UJ 

o 

o>o*^<-»r\jm 

< 

^4**—tt-tw4f* 

OC 

^mh-offiH 

00  00 

lu>z 

•    ••••• 

UJ 

•    ••••• 

•    ••••• 

am 

X-JO 

LU 

oooooo 

ac 

oooooo 

o 

oooooo 

coir\ 

*-ccz 

X 

< 

o^ 

LL 

►- 

» 

LU 

^^ 

U.O.    • 

oo 

OC 

LL      • 

OC    • 

z 

UJ 

oo 

ac  • 

»— 1 

rvjmmmm 

o 

<\j<\i<\jr\j<\j<\j 

•> 

vO  nO^O  ^o  ^o^o 

o 

0CQ- 

ooooo 

< 

oooooo 

CO 

oooooo 

o 

o     > 

LU 

o 

h- 

1  1  1  1  1  1 

o> 

»-o_i 

o 

LU  LULU  LULU 

LULU  LULU  LULU 

z 

LU  LULU  LULULU 

mno 

uza 

< 

ooooo 

UJ 

oooooo 

LU 

>}-m^co%©h- 

•  i—O 

<•— LU 

o 

ooooo 

X 

oooooo 

X 

oco*t>o*?ao 

oo 

U.ZO. 

ooooo 

h- 

oooooo 

UJ 

-«ir\cgmaDr- 

i— iLU 

oo 

UJ 

o«Njao>ro 

oooooo 

-J 

COONlMflH 

H-O 

>~0C 

X 

o*o^^iNm 

LL 

sO  sO  s£)  ^O  ^0  ^0 

UJ 

vo^rgrvjirr- 

>wo 

h-J-CL 

K- 

•  ••••• 

o 

•    ••••• 

•  ••••• 

h-«0 

h-O*. 

oooooo 

oooooo 

UJ 

oooooo 

•HJO 

>o:>ic«— 

LL 

z 

o 

oujro 

mi— 3 

c 

o 

<I 

«-*       • 

KU.O 

•— • 

o 

t-U-O 

•-«    z 

z 

»- 

COO 

CO  CO  CO  •— " 

o 

rvjmmmm 

< 

LU 

-r  >»>r>r  >*>*- 

< 

LUZLUZ 

*- 

ooooo 

K 

X 

oooooo 

_JcOO 

OUJOO 

t— 

z 

H- 

i  i  l  I  1  i 

LU  LULU  LULULU 

UJ3" 

«CO<«— 1 

•— 1 

LU  LULU  LULU 

UJ 

-J»- 

O      OK 

CO 

ooooo 

■— 1 

LL 

>rmaommco 

U.ZX 

LU      O 

a 

ooooo 

OC 

o 

mmooo>roo 

uoac 

U. COLL  .2 

Q- 

ooooo 

o 

<— iaoo<— )0DO 

o 

oixaD 

ONCO<tO 

CO 

LOsj-iriLO^-m 

OOXZ 

LU        LL 

OC 

O-C^^Nfl 

OC 

oooooo 

o 

— '■    If4^h<i    1 

z>    o 

oc>oc 

< 

•    ••••• 

< 

•    ••••• 

z 

•     •••«• 

_)  a:  oo 

LU  COLLI  CO 

_J 

oooooo 

-J 

oooooo 

*■« 

oooooo 

__><co 

COZODUJ 

r> 

o 

o 

OUJ«— 

x<ro 

o 

o 

< 

OIO 

30C3< 

z 

z 

UJ 

a:  coo. 

^v-zo 

< 

< 

oc 

96 


•«>-*o<-< 

<ooo 

K             1 

LL'LLlUJUJ 

xr^>Of^ 

H'fOOCvj 

wc\j.-<po 

a. 

f*"(MfOfVi 

<; 

(S^fO 

-1 

z    •   •   • 

n 

•-000 

0 

LO    1            1 

z 

— 

00 

•»o«— '•"-• 
<ooo 

0£ 

»-     1  1 

O 

LULU  LULU 

LL 

►-LOOLO 

—* 

► 

»CLOvO 

-J 

tfi  •  •  • 

< 

coco 

s: 

Ul    1    1 

a. 

O 

z 

1  — -1-^ 

<ooo 

<ooo 

0: 

»-     1  1 

K        1     1 

0 

LLLULUUI 

LULU  LULU 

LL 

X0CLOQ0 
»-LOLOvfi 

X&>>om 
Hf-oco 

O 

*~(\ir*-o 

•  • 

mco>* 

>±~4cr\ 

Z 

z   •   •  • 

z   •  •  • 

O 

*-ooo 

mOOO 

<— 

to  1 

to  1 

h- 

D 

-J 

—   _!,_<.— ( 

«-._ irvi«— • 

O 

<ooo 

<ooc 

to 

1-     t  1 

LULU  LULU 

►-     1  1 

LULU  LULU 

LL. 

Xf-mc 

xroinst 

O 

»-rcLO— • 

LU 

NO'ON 

LT  i-i(\JLO 

to 

a 

—ir*->fr 

CvJtMO 

t- 

> 

to    •   •    • 

to   •    •    • 

Z 

t- 

ocoo 

OOOO 

UJ 

0     1  1 

Ul    1   1 

t— 1 

O 

o 

z 

— 

— 

«»o,-t_i 

—  H^ir^ 

LL 

►- 

<ooo 

<ooo 

u. 

<i 

y~ 

t-  1  1  1 

LU 

0 

LULL'LL'LU 

LULU  LULU 

o 

<— 1 

Xvoacr- 

xtof^cr 

o 

0 

K->Or^CT* 

►—00  LOO 

z 

-Hino 

<—>rnif- 

Ot 

~- 

— '>JO>$- 

>rr*-vri^ 

UJ 

0s  •-•in 

»-tf\i—< 

»-* 

X 

Z    •    •    • 

z  •  •  • 

Ct 

LU 

—000 

•—ooo 

■=> 

O 

to  t        1 

LO  1     1     ( 

o 

Z 

u. 

>— 1 

CO 

H- 

fn-H(j 

—.t— H— ll— ! 

»H  0«— • 

z 

z 

oou 

<ouo 

<ooo 

F-4 

< 

LU 

00 

t— 

t-        1 

— ' 

•■^ 

t— LL.LULU 

LULU  LULU 

LULU  LULU 

tO 

O 

Zf^oao 

X  •— <•— »■— • 

XO<N-< 

X 

>— 

.— 

<sOi^m 

»-ir\— 1—4 

►-MOO" 

f- 

LL 

>— — <r\JLTN 

•m0>t— 'in 

wOO^CO 

•-* 

CO 

LL 

tomor*- 

•"^vO1- "n 

>J-<-H— <f\! 

JC 

* 

LU 

ki:>Tm-o 

-^•-HLn 

mojfM 

O 

0  •  •  • 

to  •  •  • 

to   •   •   • 

2: 

U 

0000 

0000 

OOOO 

a 

0     1 

0     1  1 

1— 1 

a 

to 

0 

LU 

z 

•— • 

•    »QC 

•    •£* 

•  •ac 

< 

11 

a 

iSHSX. 

uo< 

oo< 

CL 

2 

ZZoj 

ZZ-U 

ZZLU 

X 

00 

c 

<OI 

<OI 

<ax 

LU 

*£ 

LL, 

t—_|LO 

k-_JtO 

f-_ltO 

97 


THETA  =  240 
STRESSES 
ACTUAL 
COMPUTED 


THETA  =  250 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.44647E  C3 


POINT  NO.  241 

LONG. 
0.24770E  01 
0.28761E  01 


DEGREES 

TANG. 
0.45455E  C3 
0.44979E  C3 


POINT    NO.    251 

LONG. 
0.22133E    01 
0.25514E    01 


SHEAR 
0.25407E    01 
0.27506E    01 


SHEAR 
0.37883E    01 
0.39775E    01 


THETA  =  260 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45222E  03 


POINT  NO.  261 

LONG. 
0.19702E  01 
0.22239E  01 


SHEAR 
0.49382E 
0.50755E 


01 
01 


THETA  =  270  DEGREES 
STRESSES  TANG. 

ACTUAL         0.45455E  C3 
COMPUTED       0.452S9E  03 


POINT    NO.    271 

LONG. 
0.17552E    01 
0.19561E    01 


SHEAR 
0.59553E    01 
0.6C426E    01 


THETA  =  280 

STRESSES 

ACTUAL 

COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45262E 


POINT  NO.  281 

LONG. 
C3       0.15748E  01 
03       0.17848E  01 


SHEAR 
0.68088E    01 
0.6R838F    01 


THETA    =    290 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E  C3 
0.45243E  C3 


POINT  NO.  291 

LONG. 
0.14345E  01 
0.17031E  01 


SHEAR 
0.74726E 

0.75850E 


01 
01 


THETA  =  300 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45347E 


POINT    NO.     301 

LONG. 
C3  0.13385E    01 

C3  0.16666E    01 


SHEAR 
0.79267E 
0.81011E 


01 
01 


THETA  =  310 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45576E  C3 


POINT  NO.  311 

LONG. 
0.12898E  01 
0.16207E  01 


SHEAR 
0.815  73E 
0.83706E 


01 
01 


THETA  =  320 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45e36E  C3 


POINT  NO.  321 

LONG. 
0.12897E  01 
0.15346E  01 


SHEAR 
0.81573E    01 
0.83472E    01 


THETA  =  330 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E 
0.459S3E 


PCINT  NO.  331 

LONG. 
C3       0.13385E  01 
C3       0.14234E  01 


SHEAR 
0.79267E 
0.8C2  72E 


01 
01 


THETA  =  340 

STRESSES 

ACTUAL 

COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45968E  03 


PCINT  NO.  341 

LONG. 
0.14345E  01 
0.13451E  01 


SHEAR 
0.74726E 
0.74556E 


01 
01 


THETA  =  350 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  C3 
0.45770E  C3 


POINT    NO.    351 

LONG. 
0.15748E    01 
0.13746E    01 


SHEAR 
0.68088E 
0.67035E 


01 
01 
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THETA  a  240  DEGREFS  POINT  NO.  24 I 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03  0.24770E  01       0.25407E  01 

COMPUTED       0.44800E  03  0.27843E  01       0.27467E  01 

THETA  =  250  DEGREES  POINT  NO.  251 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03  0.22133E  01       0.37883E  01 

COMPUTED       0.44886E  03  C.25321E  01       0.39935E  01 

THETA  =  260  UEGPFCS  POINT  NO.  261 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03  0.19702E  01       0.49382E  01 

COMPUTED       0.45C09E  03  0.22847E  01       0.51336E  01 

THFTA  =  27C  DEGREES  POINT  NO.  271 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03  0.17552E  01       0.59553E  01 

COMPUTED       0.45157E  03  0.20507E  CI       0.61327E  01 

THETA  =  280  DEGREES  POINT  NO.  281 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455F  03  0.15748E  01       0.68088E  01 

COMPUTED       0.45314E  C3  0.18395E  <U       0.69618E  01 


THETA  =  290  DEGREFS      POINT  NO.  291 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03       0.14345E  01       0.74726E  01 

COMPUTED       0.45465E  03       0.16605E  01       0.75972E  01 

THETA  =  300  DEGREES      POINT  NO.  301 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03       0.13385E  01       0.79267E  01 

COMPUTED       0.45592E  03       0.15229E  01       0.80215E  01 

THETA  =  310  DFGREFS      POINT  NO.  311 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03       0.12898E  01       0.81573E  01 

COMPUTED       C.45684E  03       0.14343E  01       0.82239E  01 

THETA  =  320  DEGREES      POINT  NO.  321 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03       0.12897E  01       0.81573E  01 

COMPUTED       0.45728E  03       0.14007E  01       0.82001E  01 

THETA  =  330  DEGREES      POINT  NO.  331 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03       0.13385E  01       0.79267E  01 

COMPUTED       0.45722E  03       0.14253E  01       0.79524E  01 

THETA  =  34C  DEGREES      POINT  NO.  341 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455E  03       0.14345E  01      "0.74726E  01 

COMPUTED       0.45666E  03       0.15C90E  01       0.74891E  01 

THETA  =  350  DEGREES      POINT  NO.  351 

STRESSES  TANG.  LONG.  SHEAR 

ACTUAL         0.45455F  03       0.15748E  01       0.68088E  01 

COMPUTED       0.45567E  03       0.16491E  01       0.68246E  01 
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THFTA  =  24C 
STRESSES 
ACTUAL 
COMPUTED 


OEGRbfcS 

TANG. 
0.45455E 
0.45166F 


POINT  NO 


03 
03 


241 

LONG. 

0.24770E  01 

0.26709E  01 


SHEAR 
0.25407E  01 
0.26847E  01 


THETA  =  250 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
0.45199F  03 


POINT  NO.  251 

LONG. 
0.22133E  01 
0.24177E  01 


SHEAR 
0.37883E  01 
0.39286E  01 


THETA  =  260 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
0.45231E  03 


POINT  NO.  261 

LONG. 
0.19702E  01 
0.21831E  01 


SHEAR 
0.49382E  01 
0.50736E  01 


THETA  =  270 

STRESSES 

ACTUAL 

COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
0.45261F  03 


POINT  NO.  271 

LONG. 
0.17552E  01 
0.19741E  01 


SHEAR 
0.59553E  01 
0.60850E  01 


THETA  =  2H0 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
0.45288E  03 


POINT  NO.  281 

LONG. 
0.15748E  01 
0.17972E  01 


SHEAR 
0.68088E  01 
0.69321E  01 


THETA  =  290 
STRESSES 
ACTUAL 
COMPUTED 


THETA  =  300 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455F  03 
0.45311E  03 


POINT  NO.  291 

LONG. 
0.14345E  01 
0.16577E  01 


DEGREES 

TANG. 
0.45455F  03 
0.45330E  03 


POINT  NO.  301 

LONG. 
0.13385E  01 
0.15597E  01 


SHEAR 
0.74726E  01 
0.75890E  01 


SHEAR 
0.79267E  01 
0.80358E  01 


THETA  =  310 

STRESSES 

ACTUAL 

COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
0.45345E  03 


POINT  NO.  311 

LONG. 
0.12898E  01 
0.15064E  01 


SHEAR 
0.81573E 
0.82590E 


01 
01 


THtTA  =  320 
STRESSES 
ACTUAL 
COMPUTED 


OEGREES 

TANG. 
Q.45455E 
0.45353F 


POINT  NO. 


03 
03 


321 

LONG. 
12897E  01 
14994E  01 


SHEAR 
0.81573E 
0.82518E 


01 
01 


THETA  =  330 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45356F 


POINT  NO. 


03 
03 


331 

LONG. 
13385E  01 
15388E  01 


SHEAR 
0.79267E 
0.80144E 


01 
01 


THETA  =  340 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG, 
0.45455E  03 
0.45353E  03 


POINT  NO.  341 

LONG. 
0.14345E  01 
0.16234E  01 


SHEAR 
0.74726E  01 
0.75540E  01 


THETA  =  350 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E  03 
0.45345E  03 


POINT  NO.  351 

LONG. 
0.15748E  01 
0.17507E  01 


SHEAR 
0.68088E 
0.68846E 


01 
01 
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THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 
C3 


1 
LONG. 
L7552E  01 
17558E  01 


SHEAR 
0.59553E    01 
0.59553E    01 


THETA    = 

STRESSES 

ACTUAL 

COMPUTED 


10    CEGREES 

TANG. 
0.45455E 
0.45455E 


POINT    NO. 


C3 
C3 


11 

LONG. 

0.19702E    01 

0.19707E    01 


SHEAR 
0.49382E 
0.49382E 


01 
01 


THETA  - 
STRESSES 
ACTUAL 
COMPUTED 


20  CEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 

C3 


21 
LONG. 
22133E  01 
22137E  01 


SHEAR 
0.378  84E 
0.378  84E 


01 
01 


THETA    = 

STRESSES 

ACTUAL 

COMPUTED 


30    CEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 

C3 


31 
LONG. 
0.24770E  01 
0.24774E  01 


SHEAR 
0.25407E    01 
0.25407E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


40  CEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 
C3 


41 
LONG. 
0.27534E  01 
0.27538E  01 


SHEAR 
0.123 3  IE    01 
0.12331E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


50  CEGREES 

TANG. 
0.45455E 
0.45455E 


PGINT    NO. 


C3 
C3 


51 
LONG. 
0.30341E 
0.30345E 


01 
01 


SHEAR 

•0.94609E-01 

0.946  04E-01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


60  DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT    NO. 


C3 
C3 


61 

LONG. 
0.33104E 
0.33U0E 


01 
01 


SHEAR 
0.14022E    01 
0.14D22E    01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


70  DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO 


C3 
C3 


71 
LONG. 
0.35742E 
0.35748E 


01 
01 


SHEAR 
•0.26498E 
■0.26498E 


01 
01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


80  DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 
C3 


81 
LONG. 
0.38172E 
0.38180E 


01 
01 


SHEAR 

0.37997E 

■0.379  97  E 


01 
01 


THETA  = 
STRESSES 
ACTUAL 
COMPUTED 


90  DEGREES 

TANG. 
0.45455E 
0.45455E 


POINT  NO. 


C3 
C3 


91 
LONG. 
0.40322E 
0.40331E 


01 
01 


SHEAR 
0.481  68E 
0.48168E 


01 
01 


THETA  =  100 
STRESSES 
ACTUAL 
COMPUTED 


CEGREES 

TANG. 
0.45455E 
0.45455E 


POINT    NO 


C3 
C3 


101 

LONG. 

0.42126E    01 

0.42135E    01 


SHEAR 

•0.56702E 

0.56702E 


01 
01 


THETA  =  110 
STRESSES 
ACTUAL 
COMPUTED 


DEGREES 

TANG. 
O.A5A55E  C3 
0.45A55E  C3 


POINT    NO.     Ill 

LONG. 
0.43530E    01 
0.43537E    01 


SHEAR 
0.63341E    01 
0.63341E    01 
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APPENDIX  C 

DIGITAL  COMPUTER  PROGRAM  FOR  SECTION  3  AND  NUMERICAL  RESULTS 

Subroutine  MACUTO  presented  in  this  Appendix  provides  a  method  for 
applying  the  theory  developed  in  Section  3.   Data  for  simulated  situa- 
tions can  be  obtained  from  subroutine  CANTIL  presented  in  Appendix  A. 

MACUTO  has  been  programmed  following  closely  the  development  of  the 
method.   Self  explanatory  comment  statements  are  included,  which  detail 
the  type  of  data  to  be  used  and  the  procedural  sequence  of  the  method. 
Print-outs  present  the  data  processed  and  the  computed  loading. 

Use  of  the  subroutines  is  made  by  a  brief  main  program  called 
TIBISAY.   This  program  merely  dimensions  the  subroutines,  accepts  and 
prepares  the  input  data,  and  turns  control  over  to  subroutine  MACUTO. 
The  card  arrangement  of  the  source  deck  is  similar  to  that  shown  in 
Fig.  B.l. 

Some  sample  results  of  the  application  made  of  the  programs  as 
specified  in  Section  3  are  presented  following  the  card  listings  of  the 
programs. 
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Subroutine  DSIKQ  of  the  Computer  Facility  Center  of  the  Naval  Postgraduate 
School,  to  solve  linear  simultaneous  equations. 
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APPENDIX  E 

Effects  of  Changing  the  Number  of  Fourier  Coefficients  and  the 
Amount  and  Quality  of  the  Data  in  the  Application  of  the  Theory  Developed 
in  Section  2. 

The  theory  developed  in  Section  2  of  this  thesis  estimates  Fourier 
coefficients  of  the  distributions  of  axial,  tangential,  and  shear  stress- 
es assuming  that  enough  data  has  been  obtained  by  placing  a  sufficient 
number  of  strain  gage  elements  around  the  external  surface  of  a  cross 
section  of  pipe.   The  correctness  of  the  estimated  coefficients  depends 
on  their  number  and  on  the  amount  and  quality  of  the  data.   Depending  on 
the  relative  orientation  of  the  gage  elements,  they  can  pick  information 
useful  for  determining  the  coefficients  of  all  the  distributions,  or,  if 
they  are  arranged  in  certain  disadvantageous  ways,  they  may  be  able  to 
provide  information  for  none,  or  only  some  of  the  coefficients. 

The  requirement  of  inferring  a  given  number  N  of  Fourier  coeffi- 
cients, demands  having,  at  a  minimum,  a  corresponding  number  M  -  M(N) 
of  valid  data  items.   If  more  data  are  available,  there  is  a  redundance 
of  data  which  permits  more  accurate  determination  of  the  N  coefficients. 

Specific  rules  for  the  necessary  number  and  disposition  of  gage  ele- 
ments to  obtain  a  specified  number  of  correct  coefficients  are  not  avail- 
able, but  they  probably  could  be  worked  out.   Some  idea  of  the  gage  ele- 
ments requirements  can  be  obtained  by  testing  actual  cases  and  comparing 
results  with  correct  ones  obtained  from  other  methods.   The  digital  com- 
puter programs  included  in  this  thesis  can  be  used  as  they  are  to  perform 
such  tests.   The  procedure  ti  follow  Ls  to  simply  specify  in  subroutine 
ZULIA  the  number  of  gage  elements,  their  angular  positions  and  angular 
orientations  around  the  cross  section,  and  the  number  of  Fourier  co- 
efficients to  be  computed.   The  print-out  of  this  program  presents  the 
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distributions  obtained  in  this  way,  shown  adjacent  to  the  theoretically 
correct  ones,  so  that  comparisons  can  be  readily  done. 

A  preliminary  survey  of  gage  elements  arrangements  was  made  using 
that  procedure,  applied  to  data  with  no  added  random  errors.   Results 
obtained  are  presented  following  these  pages.  Many  of  these  results 
might  have  been  anticipated  from  theoretical  considerations.   The  follow- 
ing nomenclature  is  used  when  applicable: 

(a)  No.  of  coefficients 

3  to  11:   the  results  presented  were  obtained  using  from 
three  to  eleven  Fourier  coefficients  per  series.   These 
results  apply  to  that  number  of  coefficients. 
7:   results  apply  only  to  seven  coefficients. 

(b)  Results 

good:  indicates  that  the  coefficients  obtained  describe 
correctly  the  tangential,  axial  and  shear  stress  distri- 
butions. 

bad:  The  computed  coefficients  do  not  describe  either 
of  the  distributions. 

only  shear  good:   good  coefficients  were  obtained  only  for 
the  shear  stress  distribution. 
Inasmuch  as  the  data  were  without  added  errors,  the  recovery  should  be 
and  actually  was  excellent  or  no  good  at  all;  accordingly,  there  is 
nothing  subjective  in  the  judgment  "good"  or  "bad". 
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OGage  elements  located  at  12  points 


equally  spaced  30°  around  the  external 
circumference  of  the  cross  section- 


gage  elements 

arrangement  No.  of  Fourier 

at  each  point         coefficients  per  series        results 


& 


T  axis 


^ 


LxL_. 


60  °y 

:6o° 


S$.°;.. 


/; 


3  to  11  good 


3  to  11  good 


90° 
^\y\  45°  3  to  11  only  shear  good 


3  to  11  bad 


3  to  11  bad 


3  to  11  good 


3  to  11  only  shear  good 


60°  3  to  11  bad 
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II) 


Gage  elements  located  at  8  points 
equally  spaced  45°  around  the  external 
circumference  of  the  cross  section. 


gage  elements 

arrangement 
at  each  point 


No.  of  Fourier 
coefficients  per  series 


results 


45°  .45 


n£s»: 


3  to  7 


good 


same 


9  and  11 


bad 


90 


3  to  7 


good 


same 


9  and  11 


bad 


v90V 


3  to  7 


only  shear  good 


same 


9  and  11 


bad 


/*5_\_ 


3  to  11 


bad 


bi— . 


3  to  11 


bad 
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Ill) 

One  gage  element  every  10°  around  the  external 
circumference  of  the  cross  section,  totaling 
36  elements.   The  first  element  at  0°  from  the 
tangential  direction,  and  every  consecutive 
one,  incrementing  this  angle  by  10°. 

from  3  to  11  coefficients:  bad  results 


IV) 


One  gage  element  every  10°  around  the  external 
circumference  of  the  cross  section,  totaling 
36  elements.   The  first  element  at  0°  from  the 
tangential  direction,  and  every  consecutive 
one,  incrementing  this  angle  by  50°. 

from  3  to  9  coefficients:   good  results 
11  coefficients:  bad  results 
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